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This paper provides an analytical description of the transport of ultrahigh energy cosmic rays 
in an inhomogeneously magnetized intergalactic medium. This latter is modeled as a collection of 
magnetized scattering centers such as radio cocoons, magnetized galactic winds, clusters or mag- 
netized filaments of large scale structure, with negligible magnetic fields in between. Magnetic 
deflection is no longer a continuous process, it is rather dominated by scattering events. We study 
the interaction between high energy cosmic rays and the scattering agents. We then compute the 
optical depth of the Universe to cosmic ray scattering and discuss the phenomological consequences 
for various source scenarios. For typical parameters of the scattering centers, the optical depth is 
greater than unity at 5 x 10 19 eV, but the total angular deflection is smaller than unity. One im- 
portant consequence of this scenario is the possibility that the last scattering center encountered by 
a cosmic ray be mistaken with the source of this cosmic ray. In particular, we suggest that part of 
the correlation recently reported by the Pierre Auger Observatory may be affected by such delusion: 
this experiment may be observing in part the last scattering surface of ultrahigh energy cosmic rays 
rather than their source population. Since the optical depth falls rapidly with increasing energy, 
one should probe the arrival directions of the highest energy events beyond 10 20 eV on an event by 
event basis to circumvent this effect. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The problem of the origin of ultrahigh energy cosmic 
rays has generally been expressed as the conjunction of 
two questions: (i) how can particles be accelerated to 
energies in excess of 10 20 eV? (ii) why is the source not 
seen in the arrival directions of the highest energy events? 
Progress on the former question has certainly been hin- 
dered by our relative lack of knowledge on acceleration 
mechanisms and high energy processes in the most pow- 
erful astrophysical objects. Regarding the latter ques- 
tion, progress has been mostly limited by the scarcity of 
experimental data at the highest energies, at least until 
very recently. 

Indeed the first results of the Pierre Auger Observa- 
tory, which have just been published, report a significant 
correlation of the arrival directions of the highest energy 
events with a catalog of active galactic nuclei (AGN) 
closer than 75 Mpc [J, This observation certainly 
marks an important step in the search for the source of 
ultrahigh energy cosmic rays. However one should not 
overinterpret the significance of these results. In par- 
ticular, the likelihood of the reported coincidence rests 
on the comparison with isotropic arrival directions, yet 
the large scale structure is known to be highly inhomo- 
geneous at least up to 75 Mpc. Since AGN are known 
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to cluster with the large scale structure, one cannot ex- 
clude at present that the observed correlation remains 
a coincidence if the source itself clusters with the large 
scale structure @|. More will be said on these data in 
Section [TV] of the present paper. 

Furthermore, there exist other (and sometimes con- 
tradictory) claims in the literature on the existence of 
correlations of ultrahigh energy cosmic ray arrival direc- 
tions with various source catalogs the strongest be- 
ing the association with BL Lacertae objects reported in 
Refs. 0, i, i| (see also Refs. @, H). Since the existing 
data is so scarce at the highest energies, the assessment 
of the statistical significance remains a difficult task. Fi- 
nally, the reported evidence for multiplets of events tends 
to suggest that the source lies in the arrival direction of 
the events clusters. However some of these clusters show 
interacting galaxies as the sole peculiar objects on the 
line of sight @, while a more recent multiplet app ears 
correlated with interacting clusters of galaxies jlOl . fl"D | . 
Taken at face value, all these claims do not allow to draw 
a clear and consistent picture of the source of ultrahigh 
energy cosmic rays. 

It is admitted that cosmic magnetic fields must play a 
key role in this puzzle, although which role exactly is also 
a question that is still seeking for an answer. And this 
source of uncertainty is in turn related to our poor knowl- 
edge of the strength and the distribution of extragalactic 
magnetic fields (see Ref. [l2l . [l3| for detailed reviews of 
existing data). There exists a rather large body of lit- 
erature on the relation between cosmic magnetic fields 
and ultrahigh energy cosmic rays. Most studies have con- 
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structed models of extragalactic magnetic fields and then 
resorted to Monte Carlo simulations in order to quantify 
the influence of these fields on the time, energy and an- 
gular images expected in large scale detectors. One must 
however underline the analytical works of Refs. [IE HH 
on cosmic ray transport in tangled extragalactic magnetic 
fields of homogeneous power, those of Refs. [IEEE] which 
discuss the particular effect of magnetic lensing and fi- 
nally Refs. [IE [IE HE] which discuss diffusive transport 
in a magnetized supercluster. 

Earlier numerical studies have addressed the phe- 
nomenology of ultrahigh energy proton propagation in 
tan gled ma gne tic fields of homogeneous power [2l|, HE 
HE IM HE HE HE HI]. There has since been a trend 
toward more realistic magnetic field configurations. For 
instance, Refs. [IE HE HE HE HE HE] have studied the 
diffusive or non-diffusive propagation in a magnetized lo- 
cal supercluster and Ref. j35[ has brought to light the 
spectral distortions induced by the interaction of ultra- 
high energy cosmic rays with a supercluster harboring 
large scale regular magnetic fields. More recently, several 
studies have attempted to model a realistic configuration 
in which the magnetic field follows the matter density 
and then studied the transport of ultrahigh energy cos- 
mic rays in the resulting structure. In order to construct 
the magnetic field, Refs. [H HE HE HE EE EE El have 
used numerical simulations of large scale structure for- 
mation involving a passive magnetic field whose strength 
was normalized to the value measured in clusters of galax- 
ies. Refs. [IE EE El) EE] have rather reconstructed the 
extragalactic magnetic field by scaling the field strength 
to the underlying density field. 

In general, these studies have assumed the magnetic 
field to be all pervading (albeit, with a more or less pro- 
nounced degree of inhomogeneity) so that magnetic de- 
flection has been modeled as a continuous process. This 
assumption has been relaxed in Ref. [ItJ which provides 
numerical simulations of cosmic ray arrival directions af- 
ter scattering with fossils of radio-galaxy lobes. Similarly, 
Ref. [48| has mentioned the possibility of discrete cosmic 
ray interactions with localized regions of enhanced mag- 
netic fields, their discussion pointing toward clusters of 
galaxies as the main scattering agents. 

This picture in which ultrahigh energy cosmic ray 
transport occurs through random discrete events is in- 
deed more likely to be valid on distance scales up to a 
few hundreds of Mpc as a consequence of the high degree 
of clustering of matter in the Universe. For instance, even 
if the magnetic field were produced in a uniform manner 
at high redshift (see Ref. [49| for a review of models of the 
origin of large scale magnetic fields), then the present-day 
magnetic field should be highly inhomogencous, as a re- 
sult of the amplification of the magnetic field in the shear 
and compressive flows associated with the formation of 
non-linear structures [IE EE EE HE HH ( see also [H for 
a general discussion). In these simulations, voids in the 
large scale structure are essentially deprived of magnetic 
field. 



Furthermore, if one attributes the origin of the ex- 
tragalactic magnetic field to pollution by a sub-class of 
galaxies, for instance starburst galaxies [IE HE HH or 
radio-galaxies [HE HE HI] , the magnetic field configura- 
tion should resemble that of a percolating process (see 
Ref. [IE] for a clear illustration). As explained further 
below, if the filling factor of the polluted regions be- 
comes comparable to that of the filaments of large scale 
structure, the filaments themselves become the scattering 
agents as in Refs. gE EE El, HH- 

The goal of the present paper is to provide an analyt- 
ical description of ultrahigh energy cosmic ray transport 
in such an inhomogeneous medium, which is modeled by 
scattering centers embedded in an unmagnetized inter- 
galactic medium. 

These scattering centers comprise the filaments just as 
the clusters of galaxies but also all possible regions of 
locally enhanced magnetic fields, such as galactic winds, 
groups of galaxies, large scale structure shocks and fos- 
sil radio-galaxy cocoons. One motivation of the present 
work is thus to make progress toward a more realistic 
magnetic field configuration which takes into account 
those localized regions of intense magnetic activity. In 
order to do so, we first sketch a census of relevant scat- 
tering centers (Section [TTJ) then analyse their respective 
influence. 

The present work is further motivated by the fact that 
Refs. [IE EE El) H3 diverge as to the conclusions they 
draw on the influence of the extragalactic magnetic fields 
on ultrahigh energy cosmic rays, even though they try 
to construct ab initio predictions for the distribution of 
these large scale magnetic fields. This difference stems 
from the uncertainty on the origin of these magnetic 
fields, not withstanding the complexity of modeling ac- 
curately the evolution of magnetic fields in the formation 
of large scale structure. Analytical tools become useful 
in this context as they allow to parametrize the influ- 
ence of such magnetic fields on the images and spectra of 
ultrahigh energy cosmic rays. This in turn will help to 
deconvolve this effect from existing and upcoming data, 
and therefore to infer useful constraints on these mag- 
netic fields. 

In the present description, magnetic deflection is no 
longer a continuous process, but is instead dominated 
by scattering events. We thus use the notion of the op- 
tical depth of the Universe to ultrahigh energy cosmic 
ray scattering and discuss the phcnomenological conse- 
quences. In particular, we show that the optical depth 
decreases very abruptly as the energy increases, because 
the source distance scale decreases due to increasing en- 
ergy losses, and because the influence of cosmic magnetic 
fields diminishes with increasing energy. 

We argue that the energy beyond which the Universe 
becomes translucent or transparent to cosmic ray scat- 
tering may be tantalizingly close to the threshold be- 
yond which experiments search for counterparts, E ~ 
4 — 6 x 10 19 eV. This could have profound consequences 
for our interpretation of existing data. For instance, if 
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most sources lie beyond the last scattering surface, one 
could mistake the scattering centers on the last scatter- 
ing surface (such as starbursts, old radio-galaxies or giant 
shock waves) with the source of ultrahigh energy cosmic 
rays. 

The phenomenological consequences thus differ widely 
from the case of continuous deflection in an all-pervading 
medium. We thus discuss in some detail the expected ef- 
fects and their relation to current and future observations 
of cosmic ray arrival directions. 

This paper is laid out as follows. In Section [Til we 
sketch a census of possible scattering centers and their 
influence on the optical depth of the Univers to cosmic 
ray scattering. We also calculate the distance to the last 
scattering surface and compare it to the expected source 
distance scale. In Section IIII| we discuss the transport 
of cosmic rays in this strongly inhomogeneous medium 
and the expected observational consequences. We no- 
tably provide sky maps of the expected optical depth up 
to different distances for our local Universe. Finally, in 
Section [IV[ we summarize our findings and comment on 
the existing data in the framework of the present model. 
The physics of the interaction of cosmic rays with scat- 
tering centers is discussed in Appendix El 

II. THE OPTICAL DEPTH OF THE UNIVERSE 
TO HIGH ENERGY COSMIC RAY SCATTERING 

A. Scattering centers in the large scale structure 

We adopt a description in which the extragalactic mag- 
netic field is inhomogeneous. If this magnetic field orig- 
inates from a sub-class of galaxies, its configuration is 
bound to follow that of the large scale structure since 
the mixing length in the Universe is small for cosmo- 
logical standards: for typical intergalactic velocities of 
~ 300km/s, the length traveled in a Hubble time is only 
~ 4Mpc. Note that the mixing length is even less in 
filaments, for which the typical dispersion of velocities is 
of order 50 km/s. 

Obviously, at a given energy, the total optical depth 
to cosmic ray scattering is dominated by the structures 
with the largest ncr, where n represents the space den- 
sity and a the cross-section of the magnetized halo. One 
should thus focus on the radio halos of radio-galaxies, the 
magnetized winds of star forming galaxies, and on larger 
scales to clusters of galaxies and filaments as well as their 
surrounding accretion shock waves. 

1. Radio halos 

Radio halos of old radio-galaxies (deemed radio ghosts) 
have been already considered as possible sites of ultrahigh 
energy cosmic ray scattering in Ref . [47[ . This study eval- 
uates their space density as n Ig ~ 10~ 2 — KU 1 Mpc -3 , 
the radius of their magnetized halos as r rg ~ 0.5— 1 Mpc 



and their magnetic field B ig ~ 1 pG. Such quasar out- 
flows have also been examined in detail in Ref. as a 
site of magnetic pollution of the intergalactic medium, 
but their results differ from those above. These lat- 
ter authors find a much lower magnetic field strength 
B rg ~ 10 -9 G, and a substantially larger extent, r rg ~ 
1 — 5 Mpc, for a comparable space density. With respect 
to the results of Ref. [HtJ , the scattering should be dom- 
inated by the sub-population of recently formed quasars 
(at redshifts z < 4), which have r Ig ~ 2 — 4 Mpc and 
B Ig ~ 3x 10~ 9 G. The main difference between these 
calculations results from the different modeling of the 
bubble evolution. The former study assumes that the 
bubble settles in pressure equilibrium in a rather dense 
and hot intergalactic environment (with p/(p) ~ 30 and 
T ~ 10 8 K) while the latter argues that the bubble ex- 
pands until its velocity matches that of the Hubble flow 
and takes the surrounding IGM to be much colder and 
less dense (T ~ 10 4 K and p/{p) = 1). Both fix the 
magnetic strength to lie at a fraction of equipartition 
with thermal energy, although this fraction to equiparti- 
tion e B = 0.5 in Ref. [13] and e B = 0.1 in Ref. [SJ; fur- 
thermore, Ref. [57} adopt cb as the equipartition fraction 
before expansion of the bubble, assuming that the mag- 
netic field then decays with expansion. This study thus 
neglects all possible further amplification mechanisms of 
B, hence their estimate (at a given es) should be consid- 
ered as a lower limit. If one instead considers eg as the 
equipartition fraction of the magnetic field at present, the 
magnetic field strength inside the bubble can be related 
to the kinetic energy of the outflow and the size of the 
bubble as follows: 

^ = bxio-'g (-) i/2 f^M 1/2 p^y 3/; 

s 1 0.1/ V 1059 ergs/ V 1M P C / 

(1) 

This latter estimate agrees with the conclusions of 
Ref. [H| which studies the degree of magnetization of 
the IGM by radio-galaxies jets and lobes. The out- 
flow energy 10 59 ergs is an average energy for a quasar 
population [47| : it corresponds to a black hole mass 
A/bh ^ 3 x 10 7 M Q , radiating L bol ~ 3 x 10 45 ergs/s 
over 10 7 yrs (5?j]. Note however that the observational 
compilation of Ref. (6(| leads to slightly higher values 
for E Tg and B Ig . These authors have observed that the 
lobes of 70% of field radio-galaxies in their sample have a 
much higher energy content ~ 10 60 — 10 61 ergs than the 
remaining 30% in clusters (about 10 58 ergs), with a typ- 
ical volume V ~ 0.03 — 0.3 Mpc 3 and inferred minimum 
energy magnetic field strengths in the range 3 — 30 pG. If 
the magnetic field is to decay as U 2 / 3 during the subse- 
quent expansion of these bubbles, the final value for i? rg 
would be of order 0.1 pG for a typical radius r rg ~ 3 Mpc 
as above. In the following, we thus consider the possible 
range of values B rg = 1 — 10 x 10~ 8 G and typical radius 
r rg ~ 1 - 3 Mpc. 

Finally, Refs. [13, H3| estimate the space density of 
quasar outflows from the observed density of quasars at 
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high redshifts and the typical duration of the quasar 
phase (taken as 10 7 yrs). Their estimate of ~ 10~ 2 — 
1CP 1 Mpc -3 agrees with the recent determinations of the 
black hole number density at low redshifts, in particu- 
lar nf> 10 7 M Q ) ~ 2-4x lO^Mpc" 3 ^ij], although 
Ref. [621 ] reports a number density that is smaller by 
about an order of magnitude. In what follows, we thus 
consider the range n rg = 3x 10~ 3 — 3 x 10~ 2 Mpc~ 3 . 



and optimistic assumptions concerning the amplification 
of -Bgw Such amplification may have been detected in 
the outflow of M82, where a magnetic field strength as 
high as 10 pG [zO] has been reported in the first 10 kpc. 
Ref. [iH has argued that the magnetic field could be am- 
plified through the Kelvin-Helmholtz instability during 
ejection. 



2. Magnetized galactic winds 

Galactic winds have been proposed as a source of mag- 
netic pollution of the intergalactic medium by various au- 
thors, see in particular [H, H3, HE] ■ Such outflows have 
been observed in different galaxies, for instance in the 
starbursting nearby dwarf galaxy M82 with wind speed 
v ~ 2000 km/s and extension ~ 10 kpc [HI, or in mas- 
sive star forming Lyman break galaxies at high redshifts 
with wind speed v ~ 1000 km/s and extending as far 
as hundreds of kpc [64[ , maybe up to ~ 1 Mpc [65| (see 
Ref. [661 ] for a review). 

Galactic winds are also a key ingredient for theoretical 
models which attempt at explainingthe metal enrich- 
ment of the intergalactic medium [5i| [6?], [H, H^]. At 
the present time, it is not clear which galaxy type (if 
any) dominates the pollution. Starburst dwarf galax- 
ies appear more akin at producing large winds, however 
they also have a smaller gaseous content and a smaller 
energetic reservoir. In the following, we use the most 
recent simulations of Ref. ,5|| which detail the prop- 
erties of galactic winds. This study shows that the 
number of wind-blowing galaxies is relatively insensitive 
to the stellar mass of the parent galaxy in the range 
10 8 M Q < M* < 10 10 Mq as a result of the opposed 
influences of wind ram pressure and amount of infalling 
material, and that this number falls at both ends of this 
mass range. At z ~ and in this mass range, the typical 
wind radius increases slowly with galaxy mass as follows: 
r gw ~ 200 kpc for M* = 10 8 M Q , r gw ~ 800 kpc for 
M* = 1O 9 M , and r gw ~ 1 Mpc for M* = 10 10 M Q . 
Overall, the contribution n gw r gw will be dominated by 
dwarf galaxies of stellar mass M* ~ 10 9 M Q . The num- 
ber density n gw of galaxies surrounded by a wind at z = 
can be derived from the filling factor / gw of the winds; 
unfortunately, this quantity appears to depend strongly 
on the model, taking values between 2 x 10~ 2 and unity. 
The median value corresponds to / gw = 0.1 — 0.2, which 
gives a density n gw ~ / g w/Vg W — 2. - 5 x 10~ 2 Mpc -3 , 
with Vg W = (47r/3)r gw the wind volume. Note that this 
number is comparable to the number density of galaxies 
of stellar mass above 10 8 — 10 9 M Q . If the filling factor 
becomes substantially larger, the galactic winds will over- 
fill the filaments in which they reside, hence the filaments 
themselves become the scattering centers. 

Concerning the strength of B gw , Ref. [55| indicates 
that most winds have a magnetic field with B gw ~ 
10~ 8 — 10~ 7 G at z — 0, the range covering conservative 



3. Clusters of galaxies 

Clusters of galaxies are rare structures in the Uni- 
verse, n cg ~ 10 _5 /i 3 Mpc~ 3 , but they are known to 
host strong magnetic fields, with -B cg |c_~ 1 — 10 ^G in 
the innermost radius r cg | c ~ 100 kpc [HI, [7lj]. Measure- 
ments of the magnetic field in the cluster outskirts are 
rather scarce as a result of the smaller electron density 
and magnetic field strength. The minimum energy in- 
terpretation of recent synchrotron data nevertheless in- 
dicates that B cg ~ 1 ^G out to r cg ~ IMpc pf3 ] - 
Theoretical expectations tend to differ. For instance, 
Ref. [4(| shows that B varies with cluster mass, and in- 
dicates that for a massive cluster B cg ~ 1 pG within 
r cg ~ 0.2 Mpc, then falls to B cg ~ 10~ 7 G within 
r cg ~ IMpc, B cg ~ 10~ 8 G within r cg ~ 2 Mpc and 
finally B cg ~ 10~ 9 G within r cg ~ 4—5 Mpc, while 
Fig. 5 of Ref. [73[ indicates more extended magnetic fields, 
with B cg ~ 1 fiG within r cg ~ 1 Mpc, then falls to 
B cg - 10~ 7 G within r cg ~ 3 Mpc, B cg - 10" 8 G 
within r cg ~ 4 Mpc and finally B cg ~ 10 -9 G within 
r cg ~ 5 Mpc. In the following, we take these two limits 
as a range for B cg and r cg . 

Note that about half of galaxies lie outside of clusters, 
hence one can treat clusters of galaxies and the above 
field radio ghosts and field galactic winds as distinct scat- 
tering centers. 



4- Filaments and walls of large scale structure 

Filaments or walls of large scale structure are not ex- 
pected to be sources of magnetic pollution per se. How- 
ever they may be pervaded with an average magnetic 
field produced in the accretion shocks surrounding them 
or generated in and ejected by the galaxies they contain, 
provided the filling factor of the resulting magnetic pol- 
lution in the filament/ wall volume is of order unity. In 
the following, we will consider both possibilities. 

If, as before, the magnetic energy density in the fila- 
ment/wall is a fraction of the thermal energy of the 
IGM, one infers a magnetic field strength: 

(2) 

pt and Tf denoting the filament baryonic density and tem- 
perature. 
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The typical length scale of a filament is If ~ 15Mpc, 
its radius rf ~ 1 — 2 Mpc, and the typical separation 
between two filaments df ~ 25 Mpc [7J] . 

During the formation of non- linear structures, shock 
waves develop as a consequence of the infall of material 
on filaments, walls and clusters of galaxies. Numerical 
simulations indicate that the typical radius of external 
shock waves around filament it is of the order of r s h — 2— 
3 Mpc [75|, [7(| ; the typical velocity of these shock waves is 
of order w sh ~ 300 — 1000 km/s. Such shock waves have 
been proposed a site of magnetic field amplification (se e 
for instance [77|) and cosmic ray acceleration (7§.l79ll80|. 

If the magnetic field in the shock wave vicinity cor- 
responds to a fraction of equipartition with the shock 
energy density pv? hl one finds: 



10" 7 G 



eg 
0.1 



1/2 



Pcxt 



1/2 



1000 km/s 



1/2 



(3) 

Note that p cxt refers to the density of infalling mate- 
rial. The estimate e_e ~ 0.1 gives the right order of mag- 
nitude for the inferred value of magnetic field strength 
^100 fiG in young supernovae remnants assuming a typ- 
ical interstellar medium density and comparable shock 
speed HHH. 

If cosmic shock waves amplify the magnetic field up to 
the value -B sn given above, one should then expect the 
filament to be endowed with a significant fraction of _B sn 
out to the shock radius. In effect, the amount of matter 
accreted through the shock in a Hubble time in units of 
the quantity of matter contained inside the structure at 
the present time can be expressed as: 



,/a< 



Poxt V s h.H 



Pii 



0.3 



rf 



t'sh 



1000 km/s / \2Mpc 



n 



Pi 



10 (pb 



(4) 



Note that the estimate Bf given in Eq. @ agrees with 
that of -B sn to within a factor of a few (even though it 
was derived through other means). 



B. Optical depth and last scattering surface for 
cosmic ray scattering 

Depending on the strength of the magnetic field in a 
halo and its coherence length, the interaction of a par- 
ticle may either lead to diffusion inside the structure, at 
sufficiently low energy, or to a weak deflection angle, at 
higher energies. The details of the interaction between a 
particle and a magnetized structure is described in detail 
in Appendix [A"l 



1. Homogeneously distributed scattering centers 

Out of simplicity, we first assume that the scattering 
centers are distributed homogeneously in the Universe 



with a typical mean free path to interaction di, where i 
refers to the type of scattering center (e.g. magnetized 
galactic wind, radio halo, filament ...). We will discuss in 
Section III CI the influence of inhomogeneity on the con- 
clusions of the discussion that follows. For scattering 
centers of density rii and cross-section cr,, d t = (jii<Ji)~ x . 
The mean free path to interaction with any scattering 
center is written d: 



1 



(•5) 



The optical depth to ultrahigh energy cosmic ray scat- 
tering over a path length I is then defined as: 



/ 

7l 



l J2 



(6) 



To make concrete estimates, assume that one type of 
scattering center dominates, with typical interaction 
length di\ 



3.1 



100 Mpc/ V 32 Mpc 



(7) 



The above fiducial value di — 32 Mpc corresponds to 
spherical scattering centers of density rij = 10~ 2 Mpc -3 
and radius fj = 1 Mpc; however it is also a typical value 
for the interaction distance to filaments of the large scale 
structure. 

The above optical depth characterizes the number of 
scatterings along a path length I but it does not provide 
information on the angular spread of the cosmic ray im- 
age on the detector. Hence it is useful to introduce an 
effective optical depth r c g, which becomes unity when the 
path length I is such that the particle has suffered a de- 
flection of order unity. If at each scattering, the squared 
deflection is noted 8&\ , then the number of scatterings to 
achieve a deflection of order unity reads 1 / 8Q\ . The scat- 
tering length i sca tt of cosmic rays in the medium, which 
corresponds to the distance over which the deflection be- 
comes of order unity, can be written as: 



1 



We thus define the effective optical depth r e ff as: 



T e ff 



(8) 



(9) 



The angular deflection can be expressed in a simple way 
as a function of the Larmor radius r^ of the particle in 
structure i, of the magnetic field coherence length Aj of 
this structure, and of the characteristic path length fj 
through the structure, which amounts to (n/2)r s for a 
sphere of radius r s or (7r/2) 2 rf for a filament of radius 
rf (see also Appendix [SJ . Using the formula provided in 



6 



Appendix [A] [in particular Eq. (|A3|) ]. one can rewrite the 
effective optical depth as: 

( 2r2 \ \ 1 
r cff ~ I nm 1 + ■ (10) 

Obviously, one always has r e g < r. One should in- 
terpret the two optical depths as follows: r < 1 (which 
implies r e fi < 1) means that the Universe is transparent 
to cosmic ray scattering on the scale I, while r > r c g > 1 
means that the Universe is opaque over this scale, i.e. 
the accumulated angular deflection is greater than unity. 
The intermediate regime, r > 1 > r G ff is interesting; it 
corresponds to a translucent situation in which cosmic 
rays suffer one to many scatterings but the accumulated 
angular deflection remains smaller than unity. 

The phenomenology of the cosmic ray signal on the de- 
tector then depends on the typical source distance, which 
should be used for I, as well as on the characteristics of 
the scattering agents described above. Assuming recti- 
linear propagation of the particles, the source distance 
scale is of order Z max , the maximal distance that a par- 
ticle of energy E can travel without loosing its energy. 
Indeed, if the source population is continously emitting 
and homogeneous (the latter being a good approximation 
on scales beyond a few hundred Mpc), the flux F{< I) re- 
ceived from sources located within a distance / increases 
as I: 

F{< I) = nsTVuHECRi , (11) 

where n s denotes the source density and -/Vtjhecr the 
number of cosmic rays emitted by a source per unit time. 
In the case of bursting sources, one finds the same scaling 
(see Ref. (3): 

F{< I) = n s Nu HECR l . (12) 

In this equation, h s should now be understood as the rate 
of bursting sources per unit time and unit volume, and 
A^uhecr as the total number of cosmic rays emitted by 
a source. 

Hence in both cases, most of the flux comes from 
sources located at distance of order l max . In the follow- 
ing, we therefore substitute Z max for I in the expression 
of the optical depth. We will discuss apart the particu- 
lar case of rare close-by sources. One can evaluate the 
distance Z max in two ways: either as the energy loss dis- 
tance E\dE/dx\ , or as the maximal distance that a 
particle can travel, assuming it has been detected with 
energy E and the maximal energy at the source is -E max . 
In the following, we use this latter definition and assume 
-E-max = 4 x 10 20 eV. The two definitions give values that 
never differ by more than 40% however, over the energy 
range 10 17 eV -> 10 20 eV. 

If particles diffuse rather than travel rectilinearly, the 
maximum distance is instead determined by y/2Dt max , 
where D denotes the diffusion coefficient and i max = 



! max /c. This will be discussed in more detail in Sec- 
tion nrcii 

We may now plot the optical depths to scattering r and 
r e ff as functions of energy. In Fig. [H we show an exam- 
ple that ignores all scattering centers except magnetized 
galactic winds, for which we assume n gw = 10~ 2 Mpc -3 , 
r gw = 0.8 Mpc, B gw = 3 • 10" 8 G and A gw = 0.05 Mpc. 
The resulting optical depth r is shown as the dashed 
(blue) line, and the effective optical depth T cff as the 
solid (red) line. The dependence of r on E actually re- 
veals the dependence of / max on E: Z max decreases sharply 
beyond a few 10 19 eV as a consequence of pion produc- 
tion on the microwave background. The dependence of 
r e ff on E is even more pronounced, since the number of 
scatterings to achieve deflection of order unity rapidly in- 
creases with energy, roughly as E 2 beyond 10 18 eV here 
[see Eq. (TTU)) ]. The horizontal dotted line indicates an 
optical depth of order unity, while the vertical dotted 
lines indicate at which energy r G ff = 1 and r = 1 respec- 
tively, from left to right. As indicated on the figure, these 
lines delimit the energy ranges in which the Universe ap- 
pears opaque, translucent or transparent to cosmic ray 
scattering. Interestingly, for this example, the Universe 
is translucent at energies close to the threshold for pion 
production E GZK ~ 6 • 10 19 eV [H H. 

In Figure [U we show the optical depths for the various 
types of scattering centers, taken in turn, and for two sets 
of parameters defining their characteristics, as indicated 
in the caption. In principle, one should of course sum the 
different optical depths of the types of scattering centers. 
If, however, the pollution of magnetized winds and radio 
halos permeate the filaments and nothing else, one should 
of course only consider the filaments as the sole scattering 
agents. 

The two quantities Z sca tt & n d d are shown together with 
the maximal path length (or source distance scale) Z max in 
Fig. [3] for magnetized galactic winds as scattering agents, 
with the same parameters used to construct Fig. [TJ Fig- 
ure [3] illustrates in a different way the opaque, translu- 
cent or transparent nature of the Universe to cosmic ray 
scattering. 

One may also draw the analog of Fig.[2]for the distance 
to the last scattering surface d for the different types of 
scattering centers, as done in Fig. 



C. Inhomogeneity of the large scale structure - 
Analytic discussion 

The above results should be corrected for the pres- 
ence of inhomogeneity when the distances considered are 
smaller than the inhomogeneity length 100 Mpc. In a 
first approach, one may assume that all scattering cen- 
ters are clustered in the filaments of large scale structure. 
This affects transport in two ways: the typical distance 
to an interaction becomes of order d{ rather than di, but 
the typical deflection may be enhanced, as the probabil- 
ity of hitting more than one scattering center during the 



7 



1000.00 



» 10.00 
o 

si 
-t-i 

Oh 

o 1.00 

"3 
o 

£ 0.10 



0.01 
10 



opaque 


translucent 


transparent 


; \ 


\ 

\ 

: \ 
: \ 
: \ 
: \ 




Magnetized galactic winds 




\ z 

\ 

\ 

\ 

\ " 

\ : 



10 1 



10 la 
Energy [eV] 



10* 



10* 



FIG. 1: Optical depth to cosmic ray scattering by magnetized 
galactic winds, with n gw = 10 -2 Mpc -3 , B gw = 3 • 1CP 8 G, 
A gw = 50 kpc, and r gw = 0.8 Mpc. Solid line: optical depth 
T c ft to scattering by an angle of order unity, as defined in 
Eq. (|10p ; dashed line: optical depth r as defined in Eq. ©. 
In the energy range where r > r e g > 1, the Universe is 
opaque up to the energy loss distance ; in the range where 
t > 1 > r c ff , the Universe is translucent on this distance scale, 
meaning that cosmic rays suffer several to many scatterings 
but the total angular deflection remains below unity; finally, 
at energies where 1 > r > T e g, the Universe is transparent to 
cosmic ray scattering. 



interaction with a filament is itself increased. 

As the density of scattering centers in a filament be- 
comes 7ij|f = n,//f, where /f ~ 5% is the average fila- 
ment filling factor in the Universe, the mean free path to 
interaction inside a filament becomes f[di. Consequently, 
the average number of interactions iV int | f with scattering 
centers of type i during the ballistic crossing of a filament 
of radius rf is: 



intlf 



n 

ftdi 



(13) 



where ff is the characteristic path length of the parti- 
cle through the filament [see the discussion that follows 
Eq. §§§ and Appendix [A] . This formula assumes that 
the particle suffers a deflection angle much smaller than 
unity at each interaction. The particle thus exits the fila- 
ment with a total deflection and time delay (with respect 
to straight line crossing): 



= N; 



intlf' 



9 2 



2 r l 



(14) 
(15) 



In these equations, 89i and 8ti denote respectively the 
deflection angle and time delay consecutive to an inter- 
action with scattering center of type i, as discussed in 
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FIG. 2: Optical depth to cosmic ray scattering for differ- 
ent types of scattering agents and for two different sets of 
parameters in each case. Solid lines: optical depth r c ft to 
scattering by an angle of order unity, as defined in Eq. ((10} ; 
dashed lines: optical depth r as defined in Eq. (JS|. The 
vertical dotted line indicates E = 6 x 10 19 eV. Fossil ra- 
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df = 25 Mpc, Bt = 3 • 10" 9 G, A f = 300 kpc, and r f = 2 Mpc 
(lower curves); d t = 25 Mpc, B s = 3 • 10" 8 G, A f = 300 kpc, 
and rf = 2 Mpc (upper curves). 



Appendix IA1 while SOm and Sti\{ give the corresponding 
deflection angle and time delay after the crossing of a 
filament. 

In the opposite diffusive regime, in which S6f ~ 1, the 
particle follows a random walk. If the interaction length 
ffdi in the filament is much smaller than the filament 
radius rf, then the analysis of diffusive propagation in a 
filament conducted in Section [A 21 applies. The particle 
bounces on the filament and exits on a timescale rf /c at 
a distance ~ (ftdi/r{ ) 1 ^ 2 rf away from its point of first 
impact. 

Note that the filling factor of the magnetized halos in 
the filament is fi/ft, with /j ~ (4/3)ni<x,-7"j the average 
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FIG. 3: Distance to the last scattering surface by magnetized 



galactic winds, with n gw = 2-10 Mpc 
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10" 




A gw = 50 kpc, and r gw = 0.8 Mpc as in Fig. [T] Solid line: 
scattering length iscatt for a deflection of order unity, as de- 
fined in Eq. (JSJ; dashed line: distance d to the last scattering 
surface as defined in Eq. ([5]). The dotted (green) line in- 
dicates the maximal distance to the source Z max , which also 
gives the source distance scale. In the energy range where 
d < iscatt < /max, the Universe is opaque ; in the range where 
d < /max < iscatt, the Universe is translucent on the dis- 
tance scale / m ax, meaning that cosmic rays suffer several to 
many scatterings but the total angular deflection remains be- 
low unity; finally, at energies where l max < d < Z sca tt, the 
Universe is transparent to cosmic ray scattering. 



filling factor of scattering centers in the Universe. The 
filament becomes overfilled by the halos when /; > /f, 
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FIG. 4: Distance to the last scattering surface for different 
types of scattering agents and for two different sets of pa- 
rameters in each case. Solid lines: scattering length Z SC att 
for deflection of order unity, as defined in Eq. ([8}; dashed 
lines: distance d to the last scattering surface as defined in 
Eq. ([5]). The dotted (green) line indicates the source dis- 
tance scale /max- The vertical dotted line indicates the loca- 
tion of -Egzk. Fossil radio galaxies: n rg = 3 ■ 10~ 3 Mpc~ 3 , 
10" 8 G, A rg = 100 kpc, and r rg = 2 Mpc (lower 
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is satisfied, one needs not consider the multiple interac- 
tion scenario depicted above, as it suffices to consider the 
filaments themselves as the scattering centers. 

As mentioned above, the average distance to scattering 
is also modified if scattering centers cluster in filaments. 
It becomes di f : 
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Magnetized filaments of large scale structure: interseparation 
dt = 25 Mpc, Bt = 3 • 10" 9 G, A f = 300 kpc, and r t = 2 Mpc 
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and rt = 2 Mpc (upper curves). 
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for typical values of the parameters. One finds: 



1 - exp(-iV int | f ) 



(16) 



intlf 



1.3 



n 



2 Mpc 



ft 
0.05 



32 Mpc 



(17) 



as the denominator in this expression represents the 
probability of hitting a scattering center when the par- 
ticle hits a filament. The quantities di t {, SO^u and St^f 
suffice in principle to characterize the transport of the 
particle in this structured Universe and to derive the 
phenomenological consequences with respect to experi- 
mental data. To gauge the influence of the geometry, 
one should compare the above quantities to those ex- 
pected for a homogeneous scattering center distribution 



For the fiducial values used in Eq. (fT7j) . fi/f{ — 0.83, i.e. 
the halos barely overfill the filaments. This means that 
if n or m is larger than the quoted values, one must con- 
sider that the scattering centers are the filaments them- 
selves, with the average quantities df, ri and Bf discussed 
previously. Conversely, if Ti or rij is smaller, one must fol- 
low the above multiple interaction scheme. 

Finally, one can verify that on distance scales ^> df, 
the number of interactions (hence the angular deflection 
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and time delay) converge toward those obtained in the 
homogeneous case (at least for rectilinear propagation). 
In effect, the filling factor of filaments can be written 
in terms of rf and df as /f ~ (jr/2)rf/df, hence over a 
length scale d, the particle suffers N int ^d/d{ ~ d/di in- 
teractions. Qualitatively, the number of interactions per 
filament crossing compensates for the different distance 
between two zones of interaction (i.e. filaments). The 
effect of clustering of the scattering centers should thus 
be important on distance scales < 100 — 200 Mpc, since 
the distance between two filaments is of order 30 Mpc; 
beyond that distance, one can use the results derived in 
the homogeneous limit (Section III B 1[) . 

One cannot exclude a priori that an even more real- 
istic description of the hierarchical clustering of matter 
would produce a sophisticated law of probability for the 
interaction path length, leading to non-standard effects 
such as anomalous diffusion. A more realistic descrip- 
tion should also account for more complex distribution 
laws for the scattering center parameters. Monte Carlo 
simulations of particle propagation in a "realistic" scat- 
tering center distribution are best suited to address such 
issues and to provide quantitative estimates of the effect 
of inhomogeneity on the transport. 

In the following section, we describe the simulations we 
have performed in order to study the influence of a real- 
istic spatial distribution of scattering centers. In view of 
the uncertainties surrounding the origin of extragalactic 
magnetic fields and the parameters describing the scat- 
tering centers, we simply describe these latter with aver- 
age values, as discussed in Section HTl 



D. Inhomogeneity of the large scale structure - 
Numerical simulations 

We have performed our simulations using a variant of 
the numerical code described in Ref. [46| . The simula- 
tion of the dark matter density field has been produced 
by the RAMSES code (8|| , and was kindly provided to us by 
S. Colombi; its characteristics are 256 3 cells, with extent 
280 Mpc, giving a grid size 1.1 Mpc. For each simulation, 
we sample a population of scattering centers. We adopt 
two physically motivated bias models: in the first model, 
the scattering center density is proportional to the dark 
matter density field; in the second, the same proportion- 
ality applies, but we do not allow scattering centers to 
reside in regions with dark matter density p < 0.5(p). 
This latter model enhances the segregation of scattering 
centers in the large scale structure. 

Figure [5] shows an example of a scattering center dis- 
tribution in a two-dimensional slice of the simulation box 
in the first bias model. The segregation of scattering cen- 
ters in filaments of the large scale structure is apparent, 
although some tend to reside in smaller density regions 
as a result of the large volume fraction occupied by such 
regions. Out of simplicity, each scattering center is mod- 
elled as a cube of the size of a cell of the simulation; 
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FIG. 5: Distribution of scattering centers in the large scale 
structure (in white), in a model in which the density of scat- 
tering centers follows that of dark matter (density contrast 
represented in colors). The thickness of this slice is 1.1 Mpc. 
The average density is 10 -2 Mpc -3 . 



each cell in the simulation is thus occupied by zero or 
one scattering center. 

We then follow the trajectories of cosmic rays of vari- 
ous energies, using the method of Ref. |5f|, which simu- 
lates the transport of particles across cells of coherence of 
the magnetic field in both the diffusive and non-diffusive 
regime. These simulations allow to compute the various 
statistical properties of transport. A first effect brought 
to light by these simulations is the general increase in 
the length of first interaction in the inhomogeneous case, 
when compared to the homogeneous scattering center 
distribution. This increase is of order 40% for the first 
bias model, and about 60% for the second bias model. 
It does not seem to depend strongly on the scattering 
center density. 

Another significant effect is related to the source en- 
vironment. If this latter is dense, as one might expect, 
the local scattering center density is higher than aver- 
age, and therefore the cosmic ray may experience sev- 
eral interactions in the source environment in the first 
megaparsecs. Accordingly, the probability distribution 
for the first interaction departs from a simple exponen- 
tial law: it exhibits a peak in the first Mpc, then de- 
creases as an exponential. These extra interactions will 
not affect strongly the total deflection angle as seen from 
the detector, since 1 Mpc seen from 100 Mpc is sub- 
tended by an angle 0.6°. The time delay associated 
to this displacement is relatively small, being of order 
~ rS 2 /(2c) ~ 180yr(r/lMpc)(<5/0.6°) 2 (r denotes here 
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the size of the structure in which the source is embedded, 
and 5 the deflection angle associated to the displacement 
within this structure). 

This effect is apparent in Fig. [5] which shows the aver- 
age number of interactions as a function of distance, for 
different energies. The dashed lines indicate the corre- 
sponding trends for a homogeneous scattering distribu- 
tion, which go to zero when the traveled distance tends 
to zero. On the contrary, the solid lines, which corre- 
spond to the simulated inhomogeneous case, depart from 
this scaling and indicate a fixed number of interactions, of 
order 2. The exact number turns out to depend on the en- 
vironment density and has a variance of order unity. For 
a source in an environment of average density, the num- 
ber of such extra interactions in the source surroundings 
is negligible. 
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FIG. 6: Average number of interactions with scattering cen- 
ters as a function of distance traveled in a time t msx (E). 
From top to bottom, solid lines correspond to different en- 
ergies (in increasing order), as indicated in the colored ver- 
sion. Dashed lines indicate the numerical results for a ho- 
mogeneous scattering center population, while solid lines cor- 
respond to the inhomogeneous case for which the scattering 
centers are distributed according to the dark matter density. 
The dotted line indicates the analytical homogeneous result 
for non-diffusive propagation. The scattering center density 
is such that di = 32 Mpc. Each scattering center is endowed 
with a magnetic field Bi = 3 • 10~ 8 G and coherence length 
\i — 100 kpc (due to the cubic geometry of the scattering cen- 
ter, this corresponds to Bi = 2.7 x 10~ s G and A; = 100 kpc 
in a spherical cell of radius 1 Mpc) . 

Figure [6] also reveals other interesting features. In 
particular, one can see clearly that the average number 
of interactions in the inhomogeneous case converges to- 
ward that obtained in the homogeneous case on distance 



scales > 100 — 200 Mpc, as expected [see the discussion 
that follows Eq. (fTTj) ] . For the highest energies, namely 
E = 10 19 7 eV and E = 10 20 eV, there is a slight offset 
between the analytical prediction for di and the homoge- 
neous calculation; this difference is attributed to the cu- 
bic geometry of the scattering center. Particles of lower 
energies, in particular E = 10 19 eV, diffuse in the scatter- 
ing center distribution, as evidenced by the higher slope 
of the average number of interactions as a function of 
the traveled distance I. One can check in particular that 
Ni n t ~ (l/di) 2 as expected. At very large distances, this 
relation breaks down because the trajectory is cut after a 
time t mEOt ; hence less and less particles are able to travel 
beyond a distance ~ (rfmax) 1 / 2 ^ 2 . Similar features are 
observed in the second bias model. 

Finally, the same simulations can be used to compute 
the average deflection angle as a function of energy and 
traveled distance. This calculation is performed as fol- 
lows. At a predetermined distance I, one draws at ran- 
dom a certain number of "small spheres" positioned on 
the sphere of radius I around the source. These "small 
spheres" mimic the detectors located at distance I from 
the source. There must be a sufficient number of these 
"small spheres" to guarantee a sufficient signal, but not 
so many that they would overlap, in which case one would 
oversample the sphere of radius I. Each time a trajectory 
intersects one of these spheres, the angle between the 
particle incoming direction in this sphere and the source 
location is recorded. Iterating over the particles and the 
"small spheres" allows to reconstruct the probability dis- 
tribution of deflection angles. 

The result is shown in Fig. [7| for various energies, for 
the same inhomogeneous distribution of scattering cen- 
ters as above. Each cell is endowed with a magnetic field 
of strength Bi — 3 x 10~ 8 G and of coherence length 
Aj = 100 kpc. Since the cell is cubic, of size 1.1 Mpc, the 
deflection per interaction corresponds to that obtained 
for a spherical cell of radius 1 Mpc and magnetic field 
strength 5j = 2.7xl0~ 8 G. The values shown in Fig. \7\ 
have been computed at the following distances: 1000 Mpc 
for E = 10 19 eV, 600 Mpc for E = 10 19 3 eV, 400 Mpc 
for E = 10 19 7 eV and 90 Mpc for E = 10 20 eV. These 
distances are representative of / max hence of the source 
distance scale. At an energy E = 10 20 eV, the mean and 
median deflections are of order 3° and 2.6° respectively, 
while at E = 10 19 7 eV, they increase to 12° and 11.5°, 
and become larger at smaller energies. These values are 
about 30% smaller than those expected from the analyt- 
ical calculation, given in Eq. (f2"Tj) further below. This 
difference can stem from the slightly different number 
of interactions experienced by particles in the inhomoge- 
neous scattering center distribution, as compared to the 
homogeneous case (see Fig. [7]). The cubic geometry of 
scattering centers used in our simulation can also con- 
tribute to alter the values of the deflection angles. Ob- 
viously, these deflections could also be larger or smaller 
depending on the exact values of the scattering center 
characteristics, see discussion above. 
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FIG. 7: Histogram of deflection for different energies, as indi- 
cated. The values have been computed at different distances 
for the different energies: 1000 Mpc for E = 10 19 eV, 600 Mpc 
for E = 10 193 eV, 400 Mpc for E = 10 19 ' 7 eV and 90 Mpc 
for E — 10 20 eV. As before, the scattering center density is 
such that di = 32 Mpc. Each scattering center is endowed 
with a magnetic field Bi = 3 • 10~ 8 G and coherence length 
\i — 100 kpc (due to the cubic geometry of the scattering cen- 
ter, this corresponds to Bi — 2.7 x 10~ 8 G and Ai = 100 kpc 
in a spherical cell of radius 1 Mpc) . 



To summarize this discussion on the effect of inhomo- 
geneity, we note the following features: when the scat- 
tering centers correlate with the large scale structure, 
the probability law of first interaction and the number 
of interactions departs from those obtained in the homo- 
geneous case, at distances < 100 Mpc. The difference 
between the two cases depends on several factors: the 
source environment and the bias of the scattering center 
distribution with respect to the underlying dark matter 
distribution, in particular. It is found however that on 
large scales > 100 Mpc and in the weak deflection regime, 
one recovers the results of the homogeneous scattering 
center distribution discussed in Section Hi B II 

Extra interactions in the source environment, if suffi- 
ciently dense to be populated by scattering centers, may 
increase slightly the time delay with respect to straight 
line propagation but will not modify substantially the 
total deflection angle. In the diffusive regime, scattering 
occurs against filaments if the interaction length in the 
filament is smaller than the filament size, or against the 
scattering centers, if not. 



III. CONSEQUENCES FOR COSMIC RAY 
TRANSPORT 

In this section, we discuss the phenomenological con- 
sequences of the above model of cosmic ray transport 
with respect to the signatures of different source models, 
discussing in particular the absence or existence of coun- 
terparts. We will discuss in Section HVl the interpretation 
of existing data in the light of these consequences, and in 
particular the recent correlation announced by the Pierre 
Auger Observatory. 



A. Optically thin regime 

The optically thin regime, in which l max < d < Z sca tt, 
is trivial in terms of particle propagation: most parti- 
cles travel in straight line, without interacting in the in- 
tergalactic medium, hence one should expect to see the 
source directly in the arrival direction of the highest en- 
ergy events. However, in the case of gamma-ray burst 
sources, the spreading of arrival times through the inter- 
action with cosmic magnetic fields is essential to reconcile 
the gamma-ray burst rate with the rate of ultrahigh en- 
ergy cosmic ray detection [86| . In the absence of scatter- 
ing (hence time delay), such a bursting source would be 
essentially unobservable as the occurrence rate is much 
too low when compared to the lifetime of the experiment. 

Independently of the source scenario, there does not 
exist at present clear and unique evidence for counterpart 
identification, as discussed briefly in the introduction. 
Extra deflection could arise from an all-pervading in- 
tergalactic magnetic field or the Galactic magnetic field. 
The influence of an all-pervading intergalactic magnetic 
field has been discussed in previous works, see for in- 
stance Refs. [10, IH IH, HH for recent works. Note that 
our model of magnetized filaments and non-magnetized 
voids may be considered as an approximation to the more 
realistic magnetic field configurations derived in these 
studies. 

Concerning the influence of the Galactic magnetic field, 
existing models suggest that the typical deflection at the 
highest energies, say ~ 10 20 eV, are probably of the order 
of a few degrees [87], HH . Hence one would need to invoke 
the existence of an extended magnetized halo to provide 
sufficient deflection. Alternatively, one may consider a 
scenario in which most particles at the highest energies 
are heavy nuclei, which are more easily deflected. 



B. Translucent regime 

The intermediate regime, in which d < Z max < Z sca tt 
is interesting, because the typical deflection is smaller 
than unity, yet it could be sufficient to explain the lack 
of counterpart. 
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1. Transport 

Since the total deflection remains smaller than unity, 
one may describe the transport as near-ballistic with a 
non-zero time delay as measured relatively to straight 
line propagation. Furthermore, one may use in this case 
the time delay and deflection formulae obtained from ran- 
dom walk arguments in Ref. [14j ]. provided one accounts 
for the inhomogeneity of the magnetic field. In detail, 
at each scattering with scattering center i, the particle 
suffers and angular deflection 59i and exits with a delay 
Sti. The corresponding formulae for 59 i and SU are given 
in Eqs. (IX4l).(IMl). 

The total time delay 5t acquired over a path length I 
is given by the sum of the time delays acquired during 
each scattering as well as that resulting from the fact 
that the particle does not travel in a straight line from 
the source to the detector. If the particle is seen from 
the detector at a typical deflection angle 5a away from 
the source direction, then the time delay associated to 
this transverse displacement with respect to the line of 
sight is l5a 2 /(Ac) [89(. In the limit of large optical depth 
t > 1, this angle 5a 2 is written as [89l ] : 



5a 2 



-59 



(18) 



where 59 2 is the rms scattering angle per scattering event. 

On average, the particle interacts at every step of 
length d, with probability d/di of hitting a structure of 
type i. Then the total time delay and deflection acquired 
after traveling a path length / are: 
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(19) 
(20) 



To make simple estimates, consider the case in which one 
type of scattering event dominates the scattering history. 
Then the typical deflection angle reads (still assuming 
t > 1): 
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(21) 



where f; is the characteristic size of the scattering center 
[see after Eq. ([9]) and Appendix [X] . The optical depth 
to cosmic ray scattering is related to the distance and 
the geometrical characteristics of the scattering centers 
as in Eq. ([7]). This deflection may thus be non-negligible 
for typical parameters of the scattering centers discussed 
in the previous section. In all cases, the arrival direc- 
tion should point back to the last scattering center en- 
countered by the cosmic ray. Since scattering centers are 



highly magnetized regions, and as such are probably as- 
sociated with active objects such as radio-galaxies, one 
may be deceived by their presence on the line of sight, and 
interpret them as the source of ultrahigh energy cosmic 
rays. The smoking gun of such counterfeiting is the dis- 
tance scale to these objects: in this optically thick regime, 
most counterparts would be located at a distance scale 
d (which can be measured) significantly smaller than the 
expected distance scale Z max (which is known). 
The associated time delay reads: 
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lOQMpc 



5a 
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The second term on the r.h.s of Eq. (f2"0"|) indeed domi- 
nates largely over the first. It is easy to verify that the 
relation between 5a 2 (the rms angle between the line of 
sight to the source and the particle incoming velocity on 
the detector) and 59 2 (the rms velocity deflection angle 
per scattering) remains unchanged in the limit r < 1. 
Obviously, however, the solid angle of the source images 
cannot exceed that of the scattering center. 

Further effects related to the formation of angular im- 
ages are discussed in the following. 



2. Angular images 

The physics of the formation of angular images has 
been discussed in detail in Refs. [3, EH in the case 
of ultrahigh energy cosmic rays propagating in an all- 
pervading irregular magnetic field. In the model under 
consideration, differences may occur when the discrete- 
ness of scattering centers cannot be neglected. This oc- 
curs if t < 1, since r indicates the covering factor of the 
scattering centers. 

In the limit r S> 1, one may use the analysis of 
Refs. 0, EH provided one translates the quantities de- 
fined in these studies in terms of those relevant in the 
present case, such the scattering rate per interaction de- 
termined in Section[AJ One point of interest concerns the 
shape of the angular image. As discussed in Refs. [3 EH' 
the image will be centered on the source, and broadened 
by an angle 5a, if there are many uncorrelated paths 
through the scattering medium linking the source to the 
detector. In the present case, this condition remains un- 
changed in the limit t > 1, i.e. it reads 15a 3> Xi, with 
Ai the magnetic field coherence length of the scattering 
center. If ISa <§; Ai, the image will appear displaced 
from the true source position by an angle 5a, with a 
small dispersion. As discussed in Ref. [16[ , the distortion 
of the image does not modify (on average) the flux re- 
ceived from the source, in either limit considered above. 
This implies in particular that the presence of scattering 
centers does not modify the expected number of events, 
but only modifies the angular disposition of these multi- 
ple events. 

The intermediate regime 15a ~ A; is that were multi- 
ple images and magnetic lensing amplification effects may 
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become prominent (see Ref. [90| for a numerical demon- 
stration of magnetic lensing). However, as is unlikely 
to exceed a few hundreds of kpc, this intermediate regime 
is to be expected only in the limit of very small deflection: 



18a 



29 
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100 Mpc 



Sa 
1/7° 



A, 



100 kpc 



(23) 



This equation indeed suggests that typical angular im- 
ages should be broadened by Sa and centered on the 
source location. 

In the limit of small optical depth (r ~ 1) which 
becomes all the more relevant at the highest energies 
E ~ 10 20 eV, some noticeable differences can be ex- 
pected. Two questions of interest are: the shape of an- 
gular images, and the possible magnification or demagni- 
fication of images. As we argue, these effects depend on 
the hierarchy between the typical displacement ~ 18a 
in the scattering center plane (oriented perpendicular to 
the line of sight to the source), the size of the scattering 
center rj , as well as the typical distance between two scat- 
tering centers in this plane, which is given by (nj) -1 / 2 . 
Out of simplicity, we assume spherical scattering centers; 
we will argue that the conclusions remain unchanged for 
filaments. 

In order to study the limit r < 1 , it suffices to assume 
that there is only one scattering center on the line of 
sight to the source. We further assume that this scatter- 
ing structure is centered on the line of sight. The shape 
of the angular image is here as well determined by the ra- 
tio ISai/Xi. As we now argue, the flux received does not 
deviate from that expected in the absence of scattering, 
Fq = -/Vuhecr/(47tZ 2 ), provided the scattering center is 
larger than the image of the source, i.e. S ft > Sat, 
denoting by 8 Pi = Vi/l\ the typical apparent half open- 
ing angle of the scattering center with l\ the distance 
between the scattering center and the detector. If the 
opposite inequality holds (8 Pi < Sai), the flux from the 
source gets demagnified through scattering. This can be 
seen as follows. 

Each area element on the scattering structure can be 
assumed to dilute an incoming unidirectional flux into a 
beam of solid angle <5S1 ~ ^SOf (assuming small deflec- 
tion). As seen from the source, this defines a solid angle 
8Q\ a such that, if particles are emitted within 8Q\ a , they 
may be redirected toward the detector through scatter- 
ing. Then: 



AO 



Is 



sn 



(24) 



The ratio l\/l corresponds to the ratio between the half- 
opening angle of the cone of solid angle SQ\ S to 86 i. Ef- 
fects related to the finite size of the scattering center are 
considered further below. 

Now, of the flux impinging on the area element, only a 
fraction Sttd/8£l is diverted away toward the detector of 
solid angle Sfl^ = A&/1\ and area (this solid angle is 



measured relative to the scattering structure). One then 
finds that the flux received from the source is: 
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^Vuhecr <^d 
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= Fq min 



(25) 



In this equation, I2 = I — h represents the distance be- 
tween the source and the scattering center. The "min" 
function has been introduced in order to limit the angu- 
lar size of the source image to the minimum of the size 
produced by deflection and the size of the scattering cen- 
ter (which is seen through a solid angle irrf jl\ from the 
source). 

Thus, F = Fo if the solid angle Sfl\ s is smaller than the 
solid angle of the scattering structure as seen from the 
source, which amounts to 8ai < 8 Pi. This can be traced 
back to the compensation between a larger source image 
(which would lead to amplification) with the dilution of 
the signal into a beam of solid angle Sfl. 

If, on the contrary 8a, > 8 Pi, the source image is 
demagnified by the ratio F/F ~ Spf/Sa 2 , i.e. by the 
ratio of the solid angle of the scattering center to the solid 
angle that the source image would have if the scattering 
center had an infinite extent. One can generalize this 
result to the case of filamentary scattering centers, by 
noting that the flux gets demagnified by the ratio of the 
area of the scattering center to the projected area (on the 
scattering plane) of the beam of solid angle 8Cl\ s . Using 
previous fiducial values for the scattering centers, and 
assuming I2 = 1/2, one finds: 
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(26) 



However, this result considers only the influence of one 
scattering center on the line of sight. As the beam width 
exceeds the apparent size of the scattering center on the 
line of sight, one must take into account the possibil- 
ity that a fraction of the beam interacts with scatter- 
ing centers away from the line of sight. In the limit of 
small angle deflection, and still assuming Sai > 8 Pi, the 
flux received by the detector should be given by Eq. (j2l)|) 
above, multiplied by the number of scattering centers 
of the scattering plane intercepted by the beam of solid 
angle 8Cti s . We neglect the possible overlap of the pro- 
jected areas of the scattering centers, which corresponds 
to (nil)^ 1 / 2 > ri, or equivalently r < 1. This number of 
intercepted scattering structures can then be written as: 



nill\8VL\. 
0.96 



10- 2 Mpc~ 
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Hence the flux received from all intercepted scattering 
centers in the limit Sai > Spi is: 



Ft, 



N,F ~ rF n 



(28) 
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This result can be understood as follows: the number of 
intercepted scattering centers is the product of the sur- 
face density times the projected area (on the scatter- 
ing plane) of the beam of solid angle <5f2| s ; however, the 
demagnification factor is the ratio of the scattering cen- 
ter area to this latter, so that the total demagnification 
factor is the product of the surface density of scattering 
centers times the area of one scattering center, i.e. r. 
This argument remains unchanged for filamentary scat- 
tering centers. 

Equation (|28p gives the total demagnification of the 
flux from a source with one scattering structure on the 
line of sight, in the limits Soti > Sfli and (nj) -1 / 2 > r, 
(i.e. r < 1). Interestingly, the angular image is now 
decomposed into Ni distinct images of angular size S(3i 
each, of similar flux ~ Fqt /Ni, being separated from one 
another by an angle of order Sai. 

Note that, on average, there is neither magnification 
nor demagnification of the flux, as expected. Regar ding 
the limit r 3> 1, this effect has been discussed in Ref. [16[ 
in particular. Concerning the limit r < 1 discussed 
above, there are two possibilities. If Scti < Sfli, then 
as shown in Eq. (|25|) the flux is unchanged through scat- 
tering. If 6<Xi > 5 [3i , the flux of the source is demagnificd 
by t through scattering, but this occurs with probability 
~ t, which corresponds to the possibility of having one 
scattering structure on the line of sight. There is also a 
probability 1 — r of seeing the source directly (without 
scattering) together with echoes of flux tFq associated 
to scattering with structures off the line of sight. Hence 
the total flux is on average unchanged. Deviations from 
this average may occur in certain configurations, for in- 
stance through magnetic lensing, see Eq. (f23l above and 
Ref. [l6j . or in particular source scenarios, as discussed 
at the end of Section UlI B 51 further below. 



3. Experimental signatures for continuously emitting 
sources 

As far as continuously emitting sources are concerned, 
a possibly large angular deflection could prevent the de- 
tection of counterparts. Indeed, values such as di — 
30 Mpc and Bi = 10 -8 G suffice to produce a deflection of 
order 10° over a path length I = Z max at energy 4- 10 19 eV, 
which is a generic threshold energy used in the search for 
counterparts. The strong evolution of 66 with energy re- 
sults from the strong evolution of l max with E close to 
the threshold for pion production. This suggests that 
counterparts should be found at sufficiently high ener- 
gies, which of course asks for high statistics. 

Since the flux received from sources within distance 
I scales as I, one may expect to see the source in the 
arrival directions of a subset lo/l ma x of all events, lo be- 
ing defined as the distance at which the typical deflec- 
tion becomes comparable to the radius within which one 
searches for counterparts. This number (o/( mra should 
be smaller than unity, since if it were unity, it would 



mean that the total angular deflection for all sources is 
very small, hence that counterparts should have been de- 
tected. 



4- Experimental signatures for bursting sources 

Regarding bursting sources, and gamma-ray bursts in 
particular, Eq. (|22p shows that the typical time delay is 
sufficiently large to explain the lack of temporal associ- 
ation between cosmic ray arrival directions and gamma- 
ray bursts, as well as the continuous rate of detection of 
high energy cosmic rays. Recall indeed that one poten- 
tial difficulty of the gamma-ray burst scenario is to ex- 
plain the near continuous detection of cosmic rays at the 
highest energies ~ 10 20 eV, when the gamma-ray burst 
rate is only ~ 10 _3 yr _1 within the energy loss distance 
~ 100 Mpc. As noted by Waxman [86|, this difficulty may 
be overcome if the arrival time spread at of the highest 
energy events is sufficiently large, i.e. at > 10 3 yr at 
10 20 eV in particular. 

Following Ref. [l4j], we note that the magnitude of 
at /St is influenced by the number of different trajectories 
that the particle can follow from the source to the detec- 
tor. If indeed all particles follow the very same trajectory, 
<7f <C 5t, while if different particles may follow different 
trajectories, one should expect at ~ St. In the present 
model, Eq. (p?3"|) shows that the latter situation is much 
more likely, so that at /St ~ 1. Furthermore, broadening 
of the time signal at the highest energies is likely to be 
increased by stochastic pion production, which results in 
at/St ~ 1 [13. 

One may also calculate the number of gamma-ray burst 
sources which can contribute to the flux at a given energy 

e MM-- 
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This number of apparent gamma-ray bursts in the cosmic 
ray sky characterizes the amount of statistical fluctuation 
to expect around the mean flux at a given energy [l4T |. 
Using Eq. (|22|) . one obtains: 
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The magnitude of this number of apparent sources im- 
plies that the spectrum of ultrahigh energy cosmic rays 
should not reveal statistical fluctuations until energies as 
large as a few 10 20 eV, at least for these fiducial values 
that characterize the scattering centers. 
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5. Direction dependent effects 

Since the sources of protons with energies beyond the 
pion production threshold are bound to reside within 
100 — 200 Mpc, one may expect the optical depth of scat- 
tering centers to vary with the direction of observation, 
just as the density of matter. In order to discuss the influ- 
ence of such variation on existing and upcoming data, we 
have constructed sky maps of the matter concentration 
using the PSCz catalog of galaxies [92] which presently 
offers the most adequate survey for this task. 

The integrated column density of baryonic matter up 
to a distance I is shown in Fig. [H] for different maximal 
distances: / = 40, 80, 120, 160 Mpc (we adopt H = 
70km/s/Mpc). In order to correct for the incomplete- 
ness of the catalog, we have followed the prescriptions 
of Ref. [92j and smoothed the galaxy distribution with 
a variable gaussian filter, making use of the HEALPix li- 
brary (93| . The overall resolution of the maps is of order 
7°. 

These maps provide an estimate of the optical depth 
to cosmic ray scattering in the case in which the scat- 
tering centers are distributed as the galaxies. If their 
distribution is biased with respect to that of n s , for in- 
stance ni/(rii) = bi (n s / (n s )) , then the optical depth is 
expressed as the following function of N g /(N g ): 



fdlh (n g /(n g )) N g 



Jdln g /(n g ) (N g 



(31) 



and the prefactor (rii)ail = (r), see also Eq. ([7]). The 
quantity N g /(N g ) is that plotted in Fig. [8]for a distance 
I = 160 Mpc. This figure assumes no bias, in which case 
t = (t) N g /{N g ). Therefore, in order to read off r from 
Fig.[5]and the above formula, one should use I = 160 Mpc 
in the definition of (r) together with the inferred value of 
N g / (N g ) from Fig. [8] If the bias were not trivial, meaning 
bi{n g / (n g )) ^ 1, its main effect would be to increase the 
contrast of Fig. [5] Figures for particular situations can 
be provided upon demand. 

The above fiducial values for the scattering centers re- 
veal an important point: depending on the direction of 
observation, one may be in a regime of small optical 
depth r < 1 or large optical depth r > 1. This has 
several noteworthy consequences. 

First of all, the typical deflection angle becomes itself 
direction dependent. In particular, the values used in 
Eq. (f2"Tj) correspond to (r) = 3 and Sa cx r 1 / 2 . This 
simple scaling law along with Fig. [8] allow to estimate, as 
a function of the parameters characterizing the scattering 
centers, the typical deflection angle in different parts of 
the sky. 

A sky map of deflection angles had been provided pre- 
viously in Ref. [39] ] , using a constrained numerical simula- 
tion of the local Universe with an all-pervading (albeit in- 
homogeneous) magnetic field whose initial data was fixed 
at high redshift. One advantage of the present maps 
shown in Fig. [5] is to parametrize the expected deflection 



in terms of the properties of the scattering structures; in 
this sense, the above maps are more general. Ref. [44[ 
has also provided a similar map, using the PSCz galaxy 
catalog to construct the matter density field, and scaling 
the magnetic field to the matter density through the law 
B cx p 2 / 3 . The exponent 2/3 assumes isotropic com- 
pression of the magnetic field during structure formation 
and it seems that numerical simulations indicate a more 
sophisticated law, with an exponent closer to 1 (see dis- 
cussion in Rcfs. [46|, H2|). Ref. [44[ also reconstructs the 
galaxy density field on small scales by repopulating ran- 
domly the galaxy distribution using the density distribu- 
tion from the PSCz on larger scales, so that their map is 
influenced by this reconstruction on scales smaller than 
~ 7°. 

Following the discussion of Section IIII B 21 the flux of 
a source does not get demagnified nor magnified, up to 
possible magnetic lensing effects, as it crosses a region of 
scattering centers, provided the predicted apparent size 
of the source image does not exceed that of the scattering 
structure. It will however suffer demagnification in the 
opposite limit. Note that this does not contradict the 
fact that an isotropic distribution of sources will yield 
isotropic arrival directions on the detector provided that 
all arrival directions from the detector can be backtracked 
to infinity. Indeed, if a particular region of the sky is as- 
sociated with a particularly large angular deflection, the 
flux of any point source is diluted by deflection through 
the crossing of this structure; however, this deflection also 
opens a larger solid angle on the source plane, so that a 
larger number of sources can contribute, and both effects 
compensate each other. This fact has been discussed in 
particular in Ref. [94| with respect to ultrahigh energy 
cosmic ray propagation in the Galactic magnetic field. 

As mentioned in Ref. [94] , one loop hole of the above 
argument is the possible existence of so-called bottle or- 
bits, which do not connect the detector to infinity. How- 
ever one does not expect this effect to appear at the ultra- 
high energies under consideration in view of the (nearly) 
random and sporadic distribution of the scattering cen- 
ters throughout the Universe and in view of the random 
nature of the deflection suffered at each interaction. This 
assertion could be verified using dedicated numerical sim- 
ulations of particle propagation. 

Just as angular deflection, the time delay will depend 
on direction, as St cx r. Although the magnitude of the 
time delay (more precisely, of its variance) controls the 
number of bursting sources that can be seen at a given 
time, it does not influence the flux received as long as 
-/Vgrb ^> 1. Indeed, a larger St means a larger -/Vgrb (at 
a fixed value of at /St), but the flux of each gamma-ray 
bursts is decreased accordingly by the larger a t and both 
effects compensate each other exactly. However, if at a 
given energy -/Vgrb 1 in a certain region of the sky, 
one should observe a corresponding cut-off in the energy 
spectrum from this region of the sky, hence a reduced 
number of events. 

As a clear example of the above possibility, consider a 
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FIG. 8: Integrated galaxy column density as derived from the PSCz catalog of galaxies up to the maximal distances I = 40 Mpc, 
I = 80 Mpc, I = 120 Mpc and I — 160 Mpc from left to right and top to bottom (Mollweide projection). The contours give the 
column density N s in units of the mean column density {N s ) = (n g ) x 160 Mpc, with (n g ) the mean galaxy density. The grey 
mask indicates the regions of the sky that are not covered by the PSCz catalog [92] ]. 



region of the sky, of solid angle AO, in which the average 
optical depth to cosmic ray scattering r < 1. Then any 
source has a probability ~ r of having one scattering 
center on the line of sight, and therefore being seen if the 
time delay is sufficient. If there is no scattering center 
on the line of sight (with probability 1 — r), then the 
time delay is zero (in a first approximation), so that the 
probability of observing a source within AO and up to a 
distance / within the lifetime of an experiment Ai cxp ~ 
10 years is extremely small: 

P = ^AM 3 h GRB A 

~ 3 x itr 4 An " GRB Atcxp (32) 

0.1 str lO" 9 Mpc'V" 1 10yrs' V 1 

Note that 0.1 str corresponds to a region of half-opening 
angle ~ 10°. In practice, no source should be seen in this 
particular direction unless it resides in a highly magne- 
tized environment [see discussion after Section III Dl see 
also Eq. (|A6j) ]. As argued in Section Till B 21 one might 
see "echoes" of this source from scattering centers located 
away from the line of sight, provided ISoti > (nj/) -1 ' 2 . 
Even then, however, the total flux of these secondary 
images would be demagnified by r as compared to that 
expected from the source without scattering. 



In summary, the average flux expected in this solid 
angle AO is lower by a factor r than that expected from 
regions in which the optical depth is greater than unity. 

Conversely, if the source is not of the bursting type, 
one might see it directly in the arrival direction if this 
source lies in a hole of the foreground scattering center 
distribution. 



C. Opaque regime 

The opaque regime corresponds tor > T e g > 1. In this 
case, cosmic rays diffuse from the source to the detector 
as in a random billiard. 

The energy spectrum received from a given source is 
likely to be strongly modified by the presence of strongly 
magnetized scattering centers, as discussed in Ref. [35| . 
Roughly, one should observe a low-energy cut-off at an 
energy E c such that 59f < 1 for E > E c and 50f ~ 1 at 
lower energies. However, when one considers the energy 
spectrum received from an ensemble of sources, whose 
flux interacts with an ensemble of scattering centers, one 
should calculate the diffuse average flux in order to make 
contact with the measured spectrum. This average spec- 
trum should not differ from the spectrum corresponding 
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to rectilinear propagation if the diffusion theorem ap- 
plies [95| and magnetic horizon effects are unimportant, 

i.e. if the distance between two sources n s is smaller 
than the energy loss distance and the diffusion length. 
Otherwise, one should calculate the spectrum following 
the methods of Ref. [96| with the diffusion coefficient 
given below. 



1. Transport 

Assuming that the diffusion process obeys the normal 
law (r 2 ) = 2Dt, one may calculate the diffusion coeffi- 
cient D using random walk arguments. In particular, 
if one neglects the time spent in a magnetized struc- 
ture in the course of an interaction, the diffusion coef- 
ficient is related to the scattering length via the usual 
law: D — Zscattc, where the scattering length Z sca tt has 
been defined in Eq. ([5]) above. 

If the particle diffuses inside a structure during an in- 
teraction, then it actually gets trapped in this structure 
during a certain amount of time and exits backwards in a 
mirror- like fashion (see Appendix [A"|) . Consider for sim- 
plicity a single scattering agent. One may then account 
for the effect of time trapping by counting the effective 
time taken to accomplish N steps of the random walk, 
which becomes Ndi(l + 6tic/di)/c. The correction de- 
creases D by a factor (1 + dUc/di). Since the trapping 
time SU ~ Ti I c is smaller than the typical distance di be- 
tween two scattering centers, this correction is not dom- 
inant. Concerning the effect of mirroring, it suffices to 
note that it takes two interactions to achieve isotropic de- 
flection, hence this decreases the diffusion coefficient by 
another factor of 2. These two corrections thus remain 
of order unity. 

The general scaling of this diffusion coefficient with 
energy is easily grasped. At low energies (typically E < 
10 18 eV depending on the parameters characterizing the 
scattering agents), it does not depend on energy, as Z sca tt 
simply corresponds to the mean free path for scattering 
d. In the high energy regime, D oc E 2 since the number 
of scatterings to achieve a deflection of order unity scales 
in the same way. The above diffusion coefficient may be 
used to describe the propagation of particles, as done in 
Refs. ElUSSiSlHi- One may add that the influence 
of any putative all-pervading magnetic field Bigm may 
be safely neglected, even at energies of order 10 18 eV, 
as long as -Bigm S 10~ n G, since the Larmor radius 
r L ~ 100 Mpc (£/10 18 eV) (Bigm/10 -11 G)~\ 

In principle, a realistic distribution of magnetic field 
cells inside the large scale structure might induce a scat- 
tering law with a more complex profile than the standard 
exponential form adopted here, which would furthermore 
depend on time in a non-trivial way so as to account for 
the effect of trapping. The particle would then follow 
a so-called continuous time random walk with waiting 
times, the properties of which can be derived by following 
the methods developed in Refs. [9^, llOOj . It would cer- 



tainly be particularly interesting if anomalous diffusion 
laws were to occur in such magnetic field configurations. 



2. Experimental signatures for continuously emitting 
sources 

The arrival direction of high energy events will point 
back to the source only if this latter is located at a dis- 
tance closer than l sca tt- In the diffusive regime, the source 
distance scale is no longer Z max but VXscatt^max, since this 
latter gives the distance that a particle can cross before 
losing its energy. Since we assume Z max > Z sca tt, most of 
the sources are located beyond Z sca tt- 

In the steady state regime, the diffusive flux received 
from a source at distance I scales as l/(? S cattO> ncncc 
the flux received from sources within I, with / > ? S catt> 
scales as I 2 /Z sca tt- Consequently, the fraction of the flux 
that can be received from sources at distances closer than 
^ sca tt [given by Eq. flTTJ|] is roughly Z sca tt/Zmax, just as in 
the non-diffusive regime. This fraction gives the fraction 
of events behind which one can hope to detect the source. 

Note that the same delusive effect of finding a scatter- 
ing center in the arrival direction of cosmic rays occurs 
in this regime just as in the translucent regime. 

On general grounds, one expects the number of multi- 
plets to be significantly smaller in this case than for small 
deflection, since the angular size of the image is consid- 
erably broadened. However, sources within the sphere of 
large angular scattering (for which the Universe appears 
translucent) may produce images with higher multiplicity 

— 1/3 

if they exist, i.e. if n s < t sca tt- The number of events 
expected from a source at distance I can be written as: 
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(33) 



In order to derive this estimate, it suffices to express the 
flux received from this source, and to replace Auhecr i n 
this expression using Eq. (fTTjl . The parameter f cov corre- 
sponds to the sensitivity of the detector in the direction 
of the source, normalized to the average sensitivity (i.e., 
on average / cov = 1). One must emphasize that the 
above equation assumes that all sources have the same 
luminosity, which may be too restrictive. 

bmcc N m cx I /I 2 , the maximum multiplicity Ni will 



be associated to the closest source at distance 



-1/3 



Ax ~ 0.1/ cov A obs - 



-1/3. 



(34) 



To provide quantitative estimates, if n s = n Si _5 x 
10~ 5 Mpc -3 , the number of events expected from the 
closest source at energies greater than 4 x 10 19 eV is a 
fraction 7 x 10~ 3 n s l^ 3 of all observed events. This num- 



— 1/3 

ber of events becomes a fraction 0.02n s _ 5 of A Q b s above 
6 x 10 19 eV. 
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Note that the expected multiplicity is the same in this 
case than that found in the absence of magnetic fields, 
since we assume the source to be within the sphere of 
large angular scattering. 

3. Experimental signatures for bursting sources 

As far as bursting sources such as gamma-ray bursts 
are concerned, most of the above results remains un- 
changed; one simply has to replace n s iVuHECR with 
n-s-^UHECR- In the present case, the typical time spread 
corresponds to the diffusive travel time, i.e. for a source 
at distance I: 

St ~ . (35) 

Therefore the number of gamma-ray bursts sources which 
can contribute to the flux at a given energy E, at any 
time, is: 

A^grb * n s ^^- . (36) 

To make concrete estimates, at 10 20 eV, Z max ~ 
95Mpc hence A^grb ~ 20 if Z scatt = 20Mpc, assuming 
h s = 10~ 9 Mpc~ 3 yr . A*grb is larger than unity, which 
implies that one should not detect significant statistical 
fluctuation in the energy spectrum and which explains 
why one can record cosmic ray events in a near contin- 
uous manner, despite the fact that close-by gamma-ray 
bursts are such rare events. 

There will of course be an energy E c where Z sca tt = 
! max , beyond which the diffusive regime will no longer 
apply. In this case, one must use the formulae given 
in Section [ill Bl for the translucent regime. Similarly, re- 
garding sources located within the sphere of large angular 
scattering, i.e. at a distance I < i sca ttj the phenomeno- 
logical consequences are those described in Section MI Al 
if / < d, or in Section IIII Bl if d < I < Z scat t . 

IV. DISCUSSION 

A. Summary of present results 

The present work has provided an analytical descrip- 
tion of ultrahigh cosmic ray transport in highly struc- 
tured extragalactic magnetic fields. The corresponding 
configuration of the extragalactic magnetic field is that of 
a collection of scattering centers, such as halos of radio- 
galaxies or starburst galaxies, or magnetized filaments, 
with a negligible magnetic field in between. Such a con- 
figuration is generally expected in scenarios in which the 
magnetic field is produced and ejected by a sub-class of 
galaxies, or generated at the accretion shock waves of 
large scale structure. Even if the magnetic field is rather 
generated at high redshift, subsequent amplification in 



the shear and compressive flows of large scale structure 
formation tends to produce a highly structured configu- 
ration, with strong fields in the filaments of galaxies and 
weak fields in the voids 0, EH ES H3 • 

In our description, transport of cosmic rays is modeled 
as a sequence of interactions with the scattering centers, 
during which the particle acquires a non-zero deflection 
angle and time delay (with respect to straight line cross- 
ing of the magnetized region), see Appendix |A1 In Sec- 
tion [Til we have sketched a list of possible scattering cen- 
ters and their characteristics (mean free path to scatter- 
ing, magnetic field, coherence length and extent). We 
have then computed the optical depth r of the Universe 
to cosmic ray scattering as a function of energy and dis- 
tance to the source, as well as the effective optical depth 
T e ff (which is defined in such a way as to become unity 
when the total angular deflection becomes unity). As dis- 
cussed in Section [TTJ the Universe can be translucent to 
cosmic ray scattering if r > 1 > T e s, meaning that the to- 
tal deflection is smaller than unity but non-zero, opaque 
if t > T c ff > 1, or even transparent if 1 > r > r c ff. For 
typical values of the scattering centers parameters, it is 
expected that the Universe be translucent or opaque on 
the source distance scale and at energies close to the pion 
production threshold. Since this energy is that gener- 
ally used by experiments as a threshold for the search for 
counterparts, the above may have important phenomeno- 
logical consequences. 

In particular, in the translucent or opaque regime, the 
closest object lying in the cosmic ray arrival direction 
should be a scattering center. Since these scattering cen- 
ters are sites of intense magnetic activity (radio-galaxies, 
starburst galaxies, shock waves, ...), they might be mis- 
taken with the source. This peculiar feature does not 
arise in models in which magnetic deflection is a con- 
tinuous process in an all-pervading magnetic field. One 
could thus conceive an "ironic" scenario, in which cosmic 
rays are accelerated in gamma-ray bursts, but scatter 
against radio-galaxies magnetized lobes, so that one in- 
terpret these latter as the source of cosmic rays because 
they are the only active objects seen on the line of sight. 
If such counterfeiting is taking place, one should observe 
that the apparent distance scale to the source (actually 
the distance to the last scattering surface) is smaller than 
the expected distance scale to the source (as determined 
by the energy losses). This offers a simple way to test for 
the above effect. 

In the translucent regime, the source image is broad- 
ened by an angle 5a which takes values of order of 
a degree at energy 10 20 eV for the fiducial values of 
the scattering structures that we considered: interac- 
tion length di ~ 30 Mpc, extent rj ~ 1 Mpc, magnetic 
field B ~ 10~ 8 G, and coherence length A, ~ 0.1 Mpc. 
The average optical depth at distance 100 Mpc is thus 
t ~ 3 for these values. Due to the uncertainties sur- 
rounding these parameters, the deflection could however 
be larger or smaller by about an order of magnitude. 
In Section IIII B 21 we have discussed effects related to 



19 



the shape of angular images when the discreteness of the 
scattering centers is taken into account. 

The inhomogeneous distribution of matter in the local 
Universe implies that this optical depth to cosmic ray 
scattering should vary with the direction of observation. 
In Section l!!! B 5| we have provided sky maps of the inte- 
grated baryonic matter density up to different distances, 
using the PSCz catalog of galaxies. These maps allow 
to estimate the fluctuation of the optical depth in differ- 
ent directions, hence that of the deflection angle, since 
8a oc r 1 / 2 . 

In our discussion, we have taken into account the in- 
homogeneous distributions of the scattering centers, see 
Sections III CI III Dl We have shown numerically that 
on path lengths longer than ~ 200 Mpc, the effect of 
inhomogeneity is negligible, as expected for a Universe 
that is homogeneous and isotropic on these scales. The 
path length to the first interaction is generally higher by 
about 40% than in the homogeneous case if the scattering 
centers distribute according to the dark matter density. 
Since scattering centers tend to concentrate in filaments 
of large scale structure, a particle may also experience 
multiple interactions upon crossing a filament, as dis- 
cussed and quantified in Section Hi CI This explains why 
the number of interactions in the inhomogeneous case 
converges toward that of the homogeneous case on long 
path lengths. 



B. Recent data from the Pierre Auger Observatory 

In its first years of operation, the Pierre Auger Ob- 
servatory has already achieved the largest aperture (in 
km 2 ■ str • yr) [l| , and it has recently released the largest 
catalog of events above 5.7 • 10 19 eV @. In this catalog, 
20 out of 27 events originate from within 3 degrees of an 
active galactic nucleus located within 75 Mpc. 

The most straightforward interpretation is to infer that 
active galactic nuclei are the sources of ultrahigh energy 
cosmic rays. However, only one of the observed counter- 
parts is of the FR-I type (Centaurus A), all others are 
more common Seyfert galaxies. From a theoretical point 
of view, this is unexpected, since these common active 
galactic nuclei do not seem to offer the required charac- 
teristics for the acceleration to ultrahigh energies |l0l| . 
Even Centaurus A, as far as its jets are concerned, does 
not appe ar to be a likely source of ultrahigh energy cos- 
mic rays 
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Furthermore, on a pur ely experimental level, Gor- 
bunov and co-authors [103j | have recently pointed out an 
anomaly in this observed correlation. Assuming that the 
AGN seen in the arrival directions of these high energy 
events are the source of ultrahigh energy cosmic rays, 
these authors have computed the expected flux using the 
known distances to these AGN. They have observed that 
the Pierre Auger Observatory has collected zero event in 
the direction to the Virgo cluster, whereas at least six 
should be expected on the basis of the large concentra- 



tion of AGN in this direction and the small distance scale 
(assuming that the cosmic rays coming from Centaurus A 
indeed originate from this object). 

Ref. [103( 1 thus argues that this observation rules out 
the possibility that AGN are the sources of ultrahigh 
energy cosmic rays, unless the cosmic rays seen in the 
direction to Centaurus A come from further away. How- 
ever, as pointed out to us during the refereeing process, 
it could also be that the absence of AGN-like source in 
Virgo is a statistical fluctuation due to the small num- 
ber of sources in the local Universe, or that all AGN-like 
sources do not have the same cosmic ray luminosity. One 
may also ponder on the possibility that the Galactic mag- 
netic field would exhibit a particular configuration in the 
direction to Virgo (which lies toward the Galactic North 
Pole), which would prevent cosmic rays from penetrat- 
ing from this direction. Hence at present, one cannot 
exclude formally that AGN are the source of ultrahigh 
energy cosmic rays, but the data of the Pierre Auger Ob- 
servatory cannot be argued to sustend this hypothesis 
strongly either. 

Another interpretation suggests that sources of ul- 
trahigh energy cosmic rays cluster with the large scale 
structure, as AGN do, hence the observed correlation 
with AGN is a coincidence. This hypothesis deserves 
to be more carefully studied, for instance by performing 
cross-correlations of the observed arrival directions with 
galaxy c atal ogs, or by following the method introduced 
in Ref. |104| . However, assuming that the sources are 



located close to the AGN which have been seen in the ar- 
rival dir ectio ns should not resolve the flux anomaly noted 
in Ref. jl03| . also it might mitigate it somewhat. 

A third interpretation is to assume that at least part 
of the observed correlation is accidental because the scat- 
tering centers on the last scattering surface cluster with 
the large scale structure, hence with AGN. This would 
alleviate this flux anomaly, since the sources would no 
longer have to be associated with the AGN distribution. 
In particular, the events seen to arise from the Centaurus 
complex might have been deflected in its vicinity. 

As mentioned previously, this scenario can be tested 
by comparing the expected source distance scale with 
the counterpart distance scale. Interestingly, both do 
not match, as the source distance scale for particles with 
observed energy 6 x 10 19 eV is of the order of 200 Mpc, 
significantly larger than the maximum distance of 75 Mpc 
for the observed counterparts. This fact has been noted 
in Ref. Q; it remained mostly unexplained, although 
it was suggested in this work that both distance scales 
would agree if the energy scale were raised by 30%. 

More quantitatively, one can calculate the probability 
that a given event with a given observed energy origi- 
nates from a certain distance, using the fraction of the 
flux contributed by sources within a certain distance at 
a certain energy. This probability law c an be calculated 
using the techniques developed in Ref. |l05l |, then tab- 
ulated. It is then possible to calculate the probability 
of seeing 20 out 27 events from a source located within 
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75 Mpc using the events energies reported in Ref. [2j]. 
This probability is small, about 3%; the mean lies at 
15 events out of 27 coming from within 75 Mpc. If one 
restricts the set of events to those that lie outside the 
Galactic plane > 12°), with 19 out of 21 seen to cor- 
relate, the probability becomes marginal, of order 0.1% 
(the mean lies at 12 out 21 within 75 Mpc). Finally, 
if one restricts oneself to the second set of events col- 
lected after May 27 2006, and on those which lie outside 
of the Galactic plane, with 9 out of 11 seen to correlate, 
the probability becomes of order 10%, with a mean at 7 
out 11 within 75 Mpc. In this latter case, the signal is 
less significant, but the statistics is also smaller. Since 
the above estimates do not take into account the uncer- 
tainty on the energy, and since they assume continuous 
instead of stochastic energy losses, these numbers should 
be taken with caution. None theless, the above estimates 
agree with those of Ref. jl06l | , which indicate that 50% of 
protons with energy E > 6 x 10 19 eV should come from 
distances less than 100 Mpc and 90% from distances less 
than 200 Mpc. 

The above discussion suggests that, unless the energy 
scale is too low or an experimental artefact is present, the 
inferred distance scale to the source appears smaller than 
the expected source distance scale. In light of the anal- 
ysis developed in the present paper, this suggests that 
part of the correlation may actually pinpoint scattering 
centers correlating with AGN rather than the source of 
ultrahigh energy cosmic rays. Said otherwise, the Pierre 
Auger Observatory may be seeing, at least partly, the 
last scattering surface of ultrahigh energy cosmic rays, 
rather than the source population. 

In order to estimate the fraction of events that are 
likely to be contaminated by such pollution, one may 
proceed as follows. Assume first that the total deflection 
imparted to the particles with energy > 6 x 10 19 eV is 
less than the 3° radius used by the Pierre Auger Obser- 
vatory for their search. One may then calculate the the 
fraction of galaxies in the PSCz catalog up to a distance 
I = 200 Mpc, weighted appropriately, which lie within 3° 
of an AGN which is itself located closer than 75 Mpc. 
The distance I — 200 Mpc is motivated by the fact that 
90% of events with energy > 6 x 10 19 eV originate from a 
distance smaller than 200 Mpc [1061 ]. One should weigh 
each galaxy with the selection function of the PSCz cat- 
alog at the distance I of this galaxy in order to correct 
for the incompleteness of the catalog; one should also 
weigh each galaxy with a factor l/l 2 to account for flux 
dilution during propagation. In the above estimate, the 
PSCz catalog is used as a tracer of the cosmic ray source 
population, and one simply calculates the probability of 
angular coincidence with the AGN sample. The num- 
ber obtained is 0.31, which suggests that 31% of events 
above 6 x 10 19 eV could correlate with the AGN, assum- 
ing that the PSCz galaxies provide an unbiased tracer 
of the cosmic ray source population and that the mag- 
netic deflection is much smaller than the search radius of 
3°. Note that this estimate does not take into account 



the effect of the magnetic field; if one were to restrict 
the angular radius to 2° in order to account for further 
possible Galactic deflection, the above fraction would be- 
come 25%. For reference, the probability that a random 
direction on the sky falls within 3° of an AGN (located 
closer than 75 Mpc) is 0.22 (becoming 0.11 for a radius 
of 2°), which therefore gives the covering factor on the 
sky of these AGN. 

In order to account for magnetic deflection, one may 
repeat the above procedure and calculate the probabil- 
ity of coincidence to within 3° of an AGN assuming that 
the event is displaced randomly by an angle 5a from the 
location of the galaxy drawn from the PSCz catalog. Of 
course, one recovers the above result 0.31 for 5a — > 0, and 
the probability 0.22 corresponding to isotropic source dis- 
tribution for 5a ~ 1 (in practice, 5a > 45° gives a prob- 
ability 0.22). Interestingly, the fraction of contaminated 
events increases as 5a becomes of order of a few degrees: 
it equals 39% for 5a = 1°, 48% for 5a = 3°, then de- 
creases, being 45% for 5a — 5° and 43% for 5a — 7°, etc. 
If the radius of the correlation with the AGN is restricted 
to 2° to allow for further deflection in the Galactic mag- 
netic field, these numbers become 21% for 5a = 1°, 29% 
for 5a = 3° and 25% for 5a = 5°. 

The above estimates indicate that, within the assump- 
tions of the above discussion, the delusion should not 
affect all events of the Pierre Auger Observatory, but 
a significant fraction nonetheless, possibly as high as 
~ 50%. Moreover, it also indicates that intergalactic 
magnetic deflection could be larger than 3° and yet pro- 
duce a relatively significant false correlation with AGN. 
If further data from cosmic ray experiments strengthen 
the observed correlation, then the present interpretation 
would fail, unless some other effects artifically enhance 
this false correlation. 

For instance, one should point out that the above frac- 
tion of contaminated events is likely to be enhanced if 
ultrahigh energy cosmic rays originate from gamma-ray 
bursts. Indeed, as discussed in Section IIIIB5[ one ex- 
pects in this case the number of events in regions of low 
foreground density to be smaller by a factor of order t 
(t being the optical depth measured in such directions) 
when compared to that coming from regions of optical 
depth greater than unity. The main reason is that a given 
source has a probability ~ r of being located behind a 
scattering center which would provide sufficient time de- 
lay for the source to become observable wih reasonable 
probability. On the contrary, a nearby gamma-ray burst 
with no scattering center on the line of sight has a negli- 
gible probability of being seen during a time span of a few 
years as a result of the small ocurrence rate. Although 
it is difficult to give a simple estimate of the magnitude 
of this effect on the amount of false correlations, one can 
easily see that it would tend to increase this fraction by 
providing more weight to regions of high foreground den- 
sity (in which AGN are more numerous). 

In Ref. [2j , the Pierre Auger Observatory has discussed 
the evolution of the probability of null hypothesis for an 
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isotropic distribution of sources with a varying search ra- 
dius, maximum AGN redshift and minimum energy (see 
Fig. 3 of Ref. 0). The minimum probability (which indi- 
cates a maximal correlation with the AGN) corresponds 
to a search radius 3.2°. This minimum can be interpreted 
as an estimate of the amount of Galactic and intergalac- 
tic magnetic deflection if one assumes that the source ex- 
actly correlates with the AGN. Interestingly, our above 
discussion suggests that this number may be a biased 
estimate and that the intergalactic deflection could be 
slightly larger. The increase of the probability of null hy- 
pothesis at larger search radii in the Pierre Auger data 
corresponds to the fact that the covering factor of the 
search area increase rapidly with search radius, being 
already 0.50 at 6°. Concerning the redshift evolution, 
one would expect in the present model that the corre- 
lation would persist to distances as large as 200 Mpc if 
the search radius is larger than the typical intergalactic 
deflection. Unfortunately, Ref. [2| does not plot this cor- 
relation beyond 100 Mpc. It would be interesting to also 
carry out this test for different search radii. 

One should emphasize that in the present interpre- 
tation, the correlation with AGN should not persist as 
the threshold energy is decreased. Indeed, the maxi- 
mum propagation distance of particles of observed energy 
4 x 10 19 eV is of order 500 Mpc, on which scale the Uni- 
verse appears isotropic. Therefore, at these energies the 
incoming flux is increasingly isotropic, and the presence 
of scattering centers on the line of sight cannot induce 
anisotropics on an isotropic sky distribution (see discus- 
sion in Section [ill B 21 as well as the discussion on the ap- 
plication of the Liouville theorem in Ref. [94] ) . The frac- 
tion of flux contribute d by the isotropic background has 
been estimated in Ref. [l07j | in the absence of extragalac- 
tic magnetic field; it reaches 83% for E > 3 x 10 19 eV, 
and 3.6% for E > 5 x 10 19 eV. The strong rise toward 
isotropy as the threshold energy decreases is thus clear. 
This effect is present, at least qualitatively, in the data of 
the Pierre Auger Observatory (see Figure 3 of Ref. Q). 

Finally, it appears that comparing the apparent source 
distance scale with the expected one, as we have done 
above, remains the most direct and simple test of the 
present interpretation. Since there is a non-negligible 
degeneracy between the expected distance scale and the 
energy calibration, it seems mandatory to obtain a cal- 
ibration through other methods that is as accurate as 
possible. 

It also appears imperative to probe the arrival direc- 
tions on an event by event basis, focussing on the most 
energetic events. In the catalog reported in Ref. Q, 
there is only one event above 10 20 eV, whose arrival di- 
rection has a relatively small super-Galactic latitude, 
bsG — —6.5°. In the above scenario, one should ex- 
pect to find a scattering center or the source on the 
line of sight, hence it should prove useful to perform a 
deep search in this direction in the radio domain, look- 
ing for traces of synchrotron emission that would attest of 
the presence of a locally enhanced intergalactic magnetic 



field. Many more events at higher energies, as e xpected 
from future detectors such as Auger North [llOj . would 
certainly help in this regard. 

A last word should be added concerning the amount 
of magnetic deflection and the source models. In partic- 
ular, it would be interesting to examine whether (and to 
what cost) the current data could be reconciled with ul- 
trahigh energy cosmic rays being accelerated in the most 
powerful AGN, which offer stronger ground than Seyfert 
galaxies for acceleration. Such a study can only be con- 
ducted through detailed Monte-Carlo simulations which 
allow for substantial scattering angles in inhomogeneous 
magnetic fields. 



As explained in Ref. 91| , gamma-ray bursts are prob- 
ably the most elusive of possible ultrahigh energy cosmic 
ray sources, as the strongest predictions are that no coun- 
terpart should be detected, that the flux should show sig- 
nificant variations around the mean at sufficiently high 
energies (a few 10 20 eV), and that multiplets of events 
should be clustered in energy. Current data do not vio- 
late any of these predictions, but it is clear that experi- 
ments with much larger aperture at the highest energies 
will be needed to test such effects. 

It is certain that much physics and astrophysics of cos- 
mic ray sources and large scale magnetic fields remain to 
be unveiled by ongoing and future detectors. 
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APPENDIX A: PARTICLE - SCATTERING 
CENTER INTERACTION 

This section describes the interaction between a par- 
ticle and a scattering center in the large scale struc- 
ture, then computes the deflection angle and the trapping 
time in the structure before its return to non-magnetized 
voids. We consider both the cylindrical geometry, which 
is representative of an interaction with a filament, and 
the spherical geometry, which we use to model the inter- 
action with a magnetized wind or cocoon. The solution 
of the diffusion equation in a planar geometry can be 
found in Ref. 0. 

This discussion assumes that the magnetic field 
strength and the diffusion coefficient are uniform in the 
scattering center. Expectations in the more general non- 
uniform case are discussed briefly at the end of this Ap- 
pendix. 



1. Interaction with a sphere or a filament: 
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general results 

If the scattering length l sc of the particle in the scat- 
tering center is much larger than the characteristic path 
length fi through the structure, the particle is simply de- 
flected by an angle 89i and emerges after a crossing time 
t ~ fi/c. The characteristic size f j should be thought of 
as the smallest length scale of the structure, i.e. (7r/2)r s 
for a sphere or (7r/2) 2 rf for a filament. The factors of 7r/2 
account for random orientation of the incoming direction. 

The deflection angle at each interaction can be com- 
puted as follows. Consider a spherical magnetized halo 
of radius r, ; magnetic field Bi and magnetic coherence 
length A.;. The magnetic scattering length l sc of a parti- 
cle of Larmor radius in this structure determines the 
length beyond which the particle has experienced a de- 
flection of order unity. Hence, if l sc -C 7j, the particle 
undergoes diffusion in the structure so that S8i ~ 0(1). 
In the following Section, it is also shown that the distance 
traveled in the structure is very small as compared to r, , 
so that escape actually takes place close to the point of 
entry with a mirror-like deflection of order tt (to within 

±7T/2). 

If, however, Z sc 3> r il the particle is only weakly de- 
flected. In order to calculate 56?, one must specify Z sc as 
a function of and Aj . The general relationship between 
these quantities can be expressed as: 
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This equation negle cts a numerical prefactor of order 
unity (see Ref. |l02j for more details). The coefficient a 
is directly related to the level of turbulence in the struc- 
ture: 
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where SBi represents the turbulent component and Bi 
the total magnetic field. In the following, we assume 
a ~ 1, meaning full turbulence, but the calculations 
that follow may be generalized to a / 1 without diffi- 
culty. The various scenarios of magnetic pollution dis- 
cussed before do not favor the existence of significant 
coherent components of the magnetic field. 

R egard ing the exponent /?, f3 = 1 if n, 3> K [l02l 

Il08l . I lOOj ] . If, however, <C Xi, then f3 also depends 
on the shape of the turbulence spectrum. For instance, 
(3 = —2/3 for Kolmogorov turbulence, (3 = for scale 
invariant turbulence (Bohm regime). For simplicity, and 
in the absence of any knowledge of the turbulence spec- 
trum in the scattering centers, we assume (5 = 0, which 
allows to simplify the discussion. Again, it is possible to 
extend the discussion to different values of (3, albeit at 
the price of slightly more complicated expressions. 

Therefore, one finds the following deflection angle. If 
f"L <C Ai, then l sc ~ -C A, < 7"j, hence the particle 
diffuses in the structure and exits with a deflection of 



order unity. Note that the inequality Ai < simply 
states that the coherence length of the magnetic field 
cannot exceed the size of the magnetic structure. 

If tl S> Ai, then l sc ~ r^/Xi. One then must consider 
whether tl is larger or smaller than yAi^i- 111 the for- 
mer case, l sc 3> hence the particle exits with a small 
deflection angle 56f ~ fjAi/(2r£) [l^]. This numerical 
prefactor 1 /2 is valid for propagation in a turbulent mag- 
netic field; it becomes 2/3 for a randomly oriented regular 



magnetic field [14 1 



In the latter case, l sc <C r.i hence the particle exits 
with a deflection angle of order unity. 

In conclusion, the deflection angle can be written in 
the approximate form: 
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Although this form is only approximate, it interpolates 
smoothly between the two different regimes of interest 
7"l -C Vr~iXi (large deflection) and r-^ 3> \fT\Xi (small 
deflection). In the high energy (small deflection) limit, 
one finds: 
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The time delay with respect to straight line crossing 
of the magnetized structure can be calculated, following 
Refs. Ill ESS!: 
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This formula is only valid for small deflection angles; the 
corresponding time delay in the diffusive regime is dis- 
cussed further below. In the high energy limit tl 2> 
VTiXi, this gives: 
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Diffusive interaction with a filament 
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If l sc rf, the particle diffuses inside the filament 
before escaping. One can assume that the particle pen- 
etrates a length scale l sc inside the filament, and then 
enters the diffusive regime. The time-dependent diffu- 
sion equation can then be used to compute the probabil- 
ity of escape as a function of time, treating the point of 
first interaction (at depth l sc ) as an impulsive source. To 
this effect, we describe the filament as a cylinder of ra- 
dius ff and infinite extension along z and consider cylin- 
drical coordinates (r, 0,z). For simplicity, we assume a 
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spatially uniform diffusion coefficient D± in the plane 
perpendicular to z, and a spatially uniform diffusion co- 
efficient Dh in the direction along z. We also neglect en- 
ergy losses, which is justified in so far as we will show 



that the trapping time is short on the typical energy 
loss timescale. The equation for the Green's function 
g(r, 9, z, t; r , 6» , z , t Q ) reads: 



d t g - Dj-d T (rd r g) - D^djg - D {l d 2 z g = U(z - z )S(r - r )S(9 - 9 )S(t - t Q ) . (A7) 

Here, rg, 9q, Zq and to give the coordinates of the first interaction in the filament. One must also take into account 
the appropriate boundary conditions, namely that beyond radius rf, the volume is unmagnetized. In the theory of 
diffusion, such boundary conditions can be modeled by ensuring that the solution to the diffusion equation vanishes 
at a radius rf. In order to solve the diffusion equation in cylindrical coordinates under this constraint, one expands 
the angular part of the Green's function g over a basis of proper functions of the operator dg and the radial part over 
a basis of Bessel functions J m (a ms r/rf ), where a ms denotes the s— th root of J m . This guarantees that the boundary 
condition will be satisfied. The solution of Eq. (|A7[) reads: 
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The probability of having the particle inside the filament at any time t > to is then given by the volume average of g 
over the filament: 



P Ies (t;t )= dvg= V e ° s -? TT-^- ( A9 ) 
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Through the explicit decomposition of unity over the above basis of Bessel functions, one can verify that P csc (t — > 
to) = I as it should. The form of P res (i;io) tends to suggest that escape takes place on a diffusive timescale rf/D±; 
this statement is actually too naive, as shown in the following. The average residence time in the filament 5t{ is 
calculated as: 
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The factors in the sum on the r.h.s. of Eq. (|AIO|) arc much smaller than unity because the particle cannot penetrate 
further than l sc <C rf into the filament before starting to diffuse, hence ro — r{ (1 — l sc /v{). This substitution followed 
by the expansion of the Bessel functions to first order in terms of l sc /vi leads to the trapping time: 
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g =l a 0s 

which is effectively smaller than the diffusive time by a factor l sc /rf. Since D± — iZ sc c, one finally obtains: 

- E 4" = -• (A12) 

Alternatively, one could calculate the residence time by averaging Eq. (|AIO[) over the probability of first scattering 
P(ro) — exp [—(rf — r )/Z sc ] /4o but this would lead to similar results. 

We thus find that the trapping time is of order of the crossing time, an unexpected result. Since the particle 
diffuses, the linear length scale traveled in this trapping time is only / ~ Qsc/rf Y^rf <C rf. Hence the particle 
enters and exits the filament at about the same location, albeit a crossing time later. In terms of angular scattering, 
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this interaction is thus akin to mirroring, as the particle will exit in a direction separated by less than 7r/2 from the 
direction of entry in the filament. 

This law Stf ~ ff/c has been verified numerically, using Monte Carlo simulations of the interaction of a particle 
with a magnetized filament, for various coherence lengths of the magnetic field. The numerical code used has been 
described in detail in Ref. [4y|. The results are shown in Fig. [9] below, where it is seen that the average residence time 
does not depend on the coherence length of the magnetic field (which characterizes the diffusion coefficient, hence the 
scattering length), but evolves linearly with the filament radius. Numerically, one obtains Stf ~ 1.3rf/c. 



3. Diffusive interaction with a sphere 



The interaction with a sphere of radius r s is quite similar to that with a filament, although the algebra is slightly 
more cumbersome. As before, we assume that the particle penetrates a length scale Z sc before starting to diffuse in 
the magnetized sphere, and adopt appropriate boundary conditions at radius r s . The diffusion equation in spherical 
coordinates reads: 

dtg - ^d r (r 2 d r g) - -^— a d e (sm6d e g) - -y^^p = -^r~,S{r - r )S(6 - 9 )6{4> " <t>oW - to) ■ (A13) 
r l ' H sin a r 2 sin 8 r l sm 9 

Its solution is written in terms of spherical harmonics and spherical Bessel functions: 

l=+oo m—+l s— +oo ^2 D\t— to | 

g(r,9,cj),t;ro,9o,(f)o,to) = ^ ^ e ' r * Yi m (6,<p)Y lm (9o,<f)o) 

1—0 m— — l s—1 

r_o\ 2 

'■■) rlj? +1 ((3 ls ) 



xji ( Pis— ) 31 ( Pis- ) _ 3 . 2 - a , ■ (A14) 



The notation (3i s indicates the s— th zero of the spherical Bessel function jl- As before, the probability of residence 
inside the spherical structure at time t > to can be computed by integrating g over the volume: 

s=+oo 2 D\t-t \ / r \ o 
Prcs(t;to)= E e ° S ^^J0 (0OS-) R . (R s ■ (A15) 

jr[ \ r s J l3o s 3i{f3o s ) 

Here as well, P res (t — > to) = 1 as it should. Finally, the residence time can be calculated by taking the limit 
^o - * r s (l — l sc /r s ) as before and expanding to first order in l sc /r s : 

5t s = / dtP res (t;t ) ^ - Y, W = ~ ■ (A16) 

■/to C s= i Pos c 

The particle bounces on the sphere, exiting at a distance I ~ (^sc/^s) 1 ^ 2 ^ ^ r s away from its point of impact. 

I 



4. Non uniform magnetic field 

The above discussion has assumed that the magnetic 
field and the diffusion coefficient arc uniform in the scat- 
tering center. If the length scale of variation of the mag- 
netic field, l B = \ VB 2 /B 2 \- 1 is "small enough", the dis- 
cussion becomes more intricate as the scattering length 
of the cosmic ray becomes itself space dependent, and 
meaningless if it is larger than Ib- Nevertheless, one may 
expect the following to occur. 

If the scattering length as measured everywhere in the 
scattering center is larger than its size r, then the total 
deflection will remain much smaller than unity. Its value 
will be given by an average of order (tIb/v 2 ), where the 
average is to be taken on tl (through its spatial depen- 
dence via B) on the trajectory. This estimate assumes 



that the particle is deflected by 59 2 ~ (Ib/^l) 2 every Ib- 
Its corresponds to the estimate of the above discussion if 
A is replaced by I b and if B is understood as the average 
magnetic field. The crossing time will remain unchanged, 
of order r/c. 

If the scattering length is everywhere smaller than r, 
then the above results should not be modified, i.e. the 
particle will bounce on the scattering center with a trap- 
ping time of order r/c. 

Consider now the intermediate case, for instance that 
where the scattering center has a core with a magnetic 
field such that Z sc becomes smaller than the size of the 
core r c , surrounded by an envelope with B such that 
^sc ^ r- With probability ~ (r c /r) 2 , the particle may 
cross the envelope and bounce on the core; in this case the 
deflection angle is of order unity and the total crossing 
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c = 10 kpc 
c = 30 kpc 



time of order r/c. With probability ~ 1 — (r c /r) 2 , the 
particle may also cross the envelope without interacting 
with the core and suffer a deflection smaller than unity 
as calculated above; the crossing time remains the same. 
The typical deflection angle over many interactions of for 
many particles is of course given by the average of these 
two possibilities. 



5 - 



2 3 
R [Mpc] 



FIG. 9: Residence time in a magnetized filament embedded 
in a non-magnetized medium for a particle impinging on the 
filament with a scattering length l sc <C rf, as a function of 
the radius of the filament. The scattering length is a function 
of the coherence length of the magnetic field A, that corre- 
sponds to the modelling of Kolmogorov turbulence inside the 
filament, i.e. l sc oc A 2//3 , see Ref. 14611. 



[1 

[2; 
[3; 

[4; 

[5: 

[e; 
ir 

[8 

[9 

[10 
[11 



J. Abraham et al. (Pierre Auger), Science 318, 938 
(2007), arXiv:0711.2256. 

J. Abraham et al. (Pierre Auger), Astropart. Phys. 29, 

188 (2008), arXiv:0712.2843. [12[ 

T. Stanev, P. L. Biermann, J. Lloyd-Evans, J. P. [13 

Rachen, and A. A. Watson, Physical Review Letters 

75, 3056 (1995), arXiv:astro-ph/9505093. [14 

P. G. Tinyakov and I. I. Tkachev, Soviet Journal of [15 

Experimental and Theoretical Physics Letters 74, 445 

(2001), arXiv:astro-ph/0102476. 

D. S. Gorbunov, P. G. Tinyakov, I. I. Tkachev, and S. V. [16; 
Troitsky, ArXiv Astrophysics e-prints (2004), astro- 
ph/0406654. 

P. G. Tinyakov and I. I. Tkachev, PRD 69, 128301 [17 
(2004), arXiv:astro-ph/0301336. 

N. W. Evans, F. Ferrer, and S. Sarkar, PRD 69, 128302 [18 
(2004), arXiv:astro-ph/0403527. 

R. U. Abbasi, T. Abu-Zayyad, J. F. Amann, G. Arch- [19 

bold, K. Belov, J. W. Belz, S. BenZvi, D. R. Bergman, 

S. A. Blake, J. H. Boyer, et al., Astrophys. J. 636, 680 [20 

(2006). 

Y. Uchihori, M. Nagano, M. Takeda, M. Teshima, [21 
J. Lloyd-Evans, and A. A. Watson, Astroparticle 
Physics 13, 151 (2000), arXiv:astro-ph/9908193. [22[ 
G. R. Farrar, A. A. Berlind, and D. W. Hogg, Astro- 
phys. J. Lett. 642, L89 (2006), arXiv:astro-ph/0507657. 
R. U. Abbasi, T. Abu-Zayyad, J. F. Amann, G. Arch- [23 



bold, R. Atkins, J. A. Bellido, K. Belov, J. W. Belz, 
S. Y. Ben-Zvi, D. R. Bergman, et al., Astrophys. J. 623, 
164 (2005). 

P. P. Kronberg, Rep. Prog. Phys. 57, 325 (1994). 

J. P. Vallee, Fundamentals of Cosmic Physics 19, 1 

(1997). 

E. Waxman and J. Miralda-Escude, ApJ 472, 89 (1996). 
A. Achterberg, Y. A. Gallant, C. A. Norman, and D. B. 
Melrose, ArXiv Astrophysics e-prints (1999), astro- 
ph/9907060. 

D. Harari, S. Mollerach, E. Roulet, and F. Sanchez, 
Journal of High Energy Physics 3, 45 (2002), 
arXiv:astro-ph/0202362. 

D. Harari, S. Mollerach, and E. Roulet, Journal of High 
Energy Physics 7, 6 (2002), arXiv:astro-ph/0205484. 
J. Wdowczyk and A. W. Wolfendale, Nature 281, 356 
(1979). 

V. S. Berezinskii, S. I. Grigor'eva, and V. A. Dogiel, 

Astron. Astrophys. 232, 582 (1990). 

P. Blasi and A. V. Olinto, PRD 59, 023001 (1999), 

arXiv:astro-ph/9806264. 

R. Lampard, R. W. Clay, and B. R. Dawson, Astropar- 
ticle Physics 7, 213 (1997). 

G. A. Medina Tanco, E. M. de Gouveia dal Pino, and 
J. E. Horvath, Astroparticle Physics 6, 337 (1997), 
arXiv:astro-ph/9610172. 

M. Lemoine, G. Sigl, A. Olinto, and D. N. Schramm, 



26 



ApJ 486, L115 (1997). 

G. Sigl, M. Lemoine, and A. V. Olinto, PRD 56, 4470 [54 

(1997) , arXiv:astro-ph/9704204. 

R. W. Clay, S. Cook, B. R. Dawson, A. G. K. Smith, [55; 
and R. Lampard, Astroparticle Physics 9, 221 (1998). 
T. Stanev, R. Engel, A. Miicke, R. J. Protheroe, and [56 
J. P. Rachen, PRD 62, 093005 (2000), arXiv:astro- [57 ; 
ph/0003484. 

H. Yoshiguchi, S. Nagataki, S. Tsubaki, and K. Sato, [58 
Astrophys. J. 586, 1211 (2003), arXiv:astro- 
ph/0210132. [59 
O. Deligny, A. Letessier-Selvon, and E. Parizot, 
Astroparticle Physics 21, 609 (2004), arXiviastro- 
ph/0303624. [60; 
G. A. Medina Tanco, Astrophys. J. Lett. 505, L79 

(1998) , arXiv:astro-ph/9808073. [61 
G. Sigl, M. Lemoine, and P. Biermann, Astroparticle 
Physics 10, 141 (1999), arXiv:astro-ph/9806283. [62; 
M. Lemoine, G. Sigl, and P. Biermann, ArXiv Astro- 
physics e-prints (1999), astro-ph/9903124. 

Y. Ide, S. Nagataki, S. Tsubaki, H. Yoshiguchi, and 
K. Sato, Publ. Astron. Soc. Japan 53, 1153 (2001), [63 
arXiv:astro-ph/0106182. 
C. Isola, M. Lemoine, and G. Sigl, PRD 65, 023004 [64 

(2002) , arXiv:astro-ph/0104289. 

C. Isola and G. Sigl, PRD 66, 083002 (2002), 
arXiv:astro-ph/0203273. [65 
T. Stanev, D. Seckel, and R. Engel, PRD 68, 103004 

(2003) , arXiv:astro-ph/0108338. 

G. Sigl, F. Miniati, and T. A. Enfilin, PRD 68, 043002 [66 

(2003) , arXiv:astro-ph/0302388. 

G. Sigl, F. Miniati, and T. A. Enfilin, PRD 70, 043007 

(2004) , arXiv:astro-ph/0401084. 

E. Armengaud, G. Sigl, and F. Miniati, PRD 72, 043009 [67 

(2005) , arXiv:astro-ph/0412525. 

K. Dolag, D. Grasso, V. Springel, and I. Tkachev, JKAS 
37, 427 (2004). [i 
K. Dolag, D. Grasso, V. Springel, and I. Tkachev, Jour- 
nal of Cosmology and Astro-Particle Physics 1, 9 (2005), [69 
arXiv:astro-ph/0410419. 

H. Kang, S. Das, D. Ryu, and J. Cho, ArXiv e-prints 
706 (2007), 0706.2597. 

G. Sigl, PRD 75, 103001 (2007), arXiv:astro- [70 
ph/0703403. 

G. A. Medina Tanco, ArXiv Astrophysics e-prints [71 

(1997) , astro- ph/9707054. 

H. Takami, H. Yoshiguchi, and K. Sato, Astro- 
phys. J. 639, 803 (2006), arXiv:astro-ph/0506203. [72; 
H. Takami and K. Sato, ArXiv e-prints 706 (2007), 
0706.3666. 

K. Kotera and M. Lemoine (2007), arXiv:0706.1891 [73; 
[astro- ph] . 

G. Medina- Tanco and T. A. Enfilin, Astroparticle 
Physics 16, 47 (2001), arXiv:astro-ph/0011454. [74; 
V. Bcrezinsky, A. Z. Gazizov, and S. I. Grigorieva, 
ArXiv Astrophysics e-prints (2002), astro- ph/0210095. 
L. M. Widrow, Rev. Mod. Phys. 74, 775 (2002). [75 

D. Ryu, H. Kang, and P. L. Biermann, A&A 335, 19 

(1998) , arXiv:astro-ph/9803275. 

G. Sigl, F. Miniati, and T. A. Enfilin, Phys. Rev. D 70, [76 
043007 (2004). 

K. Dolag, Astronomische Nachrichten 327, 575 (2006), [77 
arXiv:astro-ph/0601484. 

P. P. Kronberg, H. Lesch, and U. Hopp, Astrophys. J. [78 



511, 56 (1999). 

G. T. Birk, H. Wiechen, H. Lesch, and P. P. Kronberg, 
Astron. Astrophys. 353, 108 (2000). 

S. Bertone, C. Vogt, and T. Enfilin, Month. Not. Roy. 
Astron. Soc. 370, 319 (2006), arXiv:astro-ph/0604462. 
M. J. Rees and G. Setti, Nature 219, 127 (1968). 
S. R. Furlanetto and A. Loeb, Astrophys. J. 556, 619 

(2001) , arXiv:astro-ph/0102076. 

Gopal-Krishna and P. J. Wiita, Astrophys. J. Lett. 560, 

L115 (2001), arXiv:astro-ph/0108117. 

S. Bertone, F. Stoehr, and S. D. M. White, Month. 

Not. Roy. Astron. Soc. 359, 1201 (2005), arXiv:astro- 

ph/0402044. 

P. P. Kronberg, Q. W. Dufton, H. Li, and S. A. Colgate, 

Astrophys. J. 560, 178 (2001), arXiv:astro-ph/0106281. 

L. Ferrarese and H. Ford, Space Science Reviews 116, 

523 (2005), arXiv:astro-ph/0411247. 

T. R. Lauer, S. M. Faber, D. Richstone, K. Gebhardt, 

S. Tremaine, M. Postman, A. Dressier, M. C. Aller, 

A. V. Filippenko, R. Green, et al., Astrophys. J. 662, 

808 (2007), arXiv:astro-ph/0606739. 

E. R. Seaquist and N. Odegard, Astrophys. J. 369, 320 

(1991). 

M. Pettini, S. A. Rix, C. C. Steidel, K. L. Adelberger, 
M. P. Hunt, and A. E. Shapley, Astrophys. J. 569, 742 

(2002) , arXiv:astro-ph/0 110637. 

K. L. Adelberger, C. C. Steidel, A. E. Shapley, and 
M. Pettini, Astrophys. J. 584, 45 (2003), arXiv:astro- 
ph/0210314. 

T. M. Heckman, in Gas and Galaxy Evolution, edited 
by J. E. Hibbard, M. Rupen, and J. H. van Gorkom 
(2001), vol. 240 of Astronomical Society of the Pacific 
Conference Series, pp. 345 — h 

A. Aguirre, L. Hernquist, J. Schaye, D. H. Weinberg, 
N. Katz, and J. Gardner, Astrophys. J. 560, 599 (2001), 
arXiv:astro-ph/0006345. 

R. Cen, K. Nagamine, and J. P. Ostriker, Astrophys. J. 
635, 86 (2005), arXiv:astro-ph/0407143. 

E. Scannapieco, C. Pichon, B. Aracil, P. Petitjean, 
R. J. Thacker, D. Pogosyan, J. Bergeron, and H. M. P. 
Couchman, MNRAS 365, 615 (2006), arXiv:astro- 
ph/0503001. 

H. -P. Reuter, U. Klein, H. Lesch, R. Wielebinski, and 
P. P. Kronberg, Astron. Astrophys. 256, 10 (1992). 

T. E. Clarke, P. P. Kronberg, and H. Bohringer, 
Astrophys. J. Lett. 547, Llll (2001), arXiv:astro- 
ph/0011281. 

F. Govoni, M. Murgia, L. Feretti, G. Giovannini, 
K. Dolag, and G. B. Taylor, Astron. Astrophys. 460, 
425 (2006), arXiv:astro-ph/0608433. 

M. Briiggen, M. Ruszkowski, A. Simionescu, M. Hoeft, 
and C. Dalla Vecchia, Astrophys. J. Lett. 631, L21 
(2005), arXiv:astro-ph/0508231. 

A. G. Doroshkevich, D. L. Tucker, R. Fong, V. Tur- 
chaninov, and H. Lin, Month. Not. Roy. Astron. Soc. 
322, 369 (2001). 

F. Miniati, D. Ryu, H. Kang, T. W. Jones, R. Cen, 
and J. P. Ostriker, Astrophys. J. 542, 608 (2000), 
arXiv:astro-ph/0005444. 

H. Kang, D. Ryu, R. Cen, and D. Song, ApJ 620, 21 
(2005), arXiv:astro-ph/0410477. 

R. M. Kulsrud, R. Cen, J. P. Ostriker, and D. Ryu, ApJ 

480, 481 (1997), arXiv:astro-ph/9607141. 

A. Loeb and E. Waxman, Nature 405, 156 (2000), 



27 



arXiv:astro-ph/0003447. 
[79] F. Miniati, MNRAS 337, 199 (2002), arXiv:astro- 
ph/0203014. 

[80] U. Keshet, E. Waxman, A. Loeb, V. Springel, and 
L. Hernquist, ApJ 585, 128 (2003), arXiv:astro- 
ph/0202318. 

[81] J. Vink and J. M. Laming, ApJ 584, 758 (2003), 

arXiv:astro-ph/0210669. 
[82] E. G. Berezhko, L. T. Ksenofontov, and H. J. Volk, AA 

412, Lll (2003), arXiv:astro-ph/0310862. 
[83] K. Greisen, Phys. Rev. Lett. 16, 748 (1966). 
[84] G. T. Zatsepin and V. A. Kuzmin, JETP Lett. 4, 78 

(1966). 

[85] R. Teyssier, A&A 385, 337 (2002). 
[86] E. Waxman, Physical Review Letters 75, 386 (1995), 
arXiv:astro-ph/9505082. 

[87] J. Alvarez-Muniz, R. Engel, and T. Stanev, Astro- 
phys. J. 572, 185 (2002), arXiv:astro-ph/0112227. 

[88] H. Takami and K. Sato, e-prints (2007), 
arXiv:0711.2386. 

[89] C. Alcock and S. Hatchett, Astrophys. J. 222, 456 
(1978). 

[90] G. Sigl, M. Lemoine, and P. Biermann, Astropart. Phys. 
10, 141 (1999), astro-ph/9806283. 

[91] E. Waxman, in Physics and Astrophysics of Ultra-High- 
Energy Cosmic Rays, edited by M. Lemoine and G. Sigl 
(2001), vol. 576 of Lecture Notes in Physics, Berlin 
Springer Verlag, pp. 122 — K 

[92] W. Saunders, W. J. Sutherland, S. J. Maddox, O. Kee- 
ble, S. J. Oliver, M. Rowan-Robinson, R. G. McMa- 
hon, G. P. Efstathiou, H. Tadros, S. D. M. White, 
et al., Month. Not. Roy. Astron. Soc. 317, 55 (2000), 
arXiv :astro-ph /000 1 1 1 7. 

[93] K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, 



F. K. Hansen, M. Reinecke, and M. Bartelmann, Astro- 
phys. J. 622, 759 (2005), arXiv:astro-ph/0409513. 

[94] D. Harari, S. Mollerach, and E. Roulet, Journal of High 
Energy Physics 8, 22 (1999), arXiv:astro-ph/9906309. 

[95] R. Aloisio and V. Berezinsky, ApJ 612, 900 (2004). 

[96] V. Berezinsky and A. Gazizov, ApJ 643, 8 (2006). 

[97] M. Lemoine, Phys. Rev. D 71, 083007 (2005). 

[98] R. Aloisio and V. Berezinsky, ApJ 625, 249 (2005). 

[99] J.-P. Bouchaud and A. Georges, Phys. Rep. 195, 127 
(1990). 

[100] R. C. Ball, S. Havlin, and G. H. Weiss, Journal of 
Physics A Mathematical General 20, 4055 (1987). 

[101] C. A. Norman, D. B. Melrose, and A. Achterberg, As- 
trophys. J. 454, 60 (1995). 

[102] F. Casse, M. Lemoine, and G. Pelletier, Phys. Rev. D 
65, 023002 (2002). 

[103] D. Gorbunov, P. Tinyakov, I. Tkachev, and S. Troitsky 
(2007), arXiv:071 1.4060 [astro- ph]. 

[104] E. Waxman, K. B. Fisher, and T. Piran, Astro- 
phys. J. 483, 1 (1997), arXiv:astro-ph/9604005. 

[105] V. Berezinsky, A. Gazizov, and S. Grigorieva, Phys. 
Rev. D 74, 043005 (2006), arXiv:hep-ph/0204357. 

[106] D. Harari, S. Mollerach, and E. Roulet, Journal of 
Cosmology and Astro-Particle Physics 11, 12 (2006), 
arXiv:astro-ph/0609294. 

[107] A. Cuoco, R. D'Abrusco, G. Longo, G. Miele, and 
P. D. Serpico, Journal of Cosmology and Astro-Particle 
Physics 1, 9 (2006), arXiv:astro-ph/0510765. 

[108] J. Giacalone and J. R. Jokipii, Astrophys. J. 520, 204 
(1999). 

[109] J. Candia and E. Roulet, JCA P 0410, 007 (2004). 
[110] http://www.augernorth.org 



